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Abstract 

A  ricochet  model  for  long  rods  developed  by  the  author  has  been  revisited  in  search  of  solution  efficiencies.  By  replacing 
a  moment-of-momentum  governing  equation  in  the  original  model  with  a  simpler  minimization  constraint,  a  solution 
technique  has  been  developed  that  avoids  the  complexity  of  an  iterative  solution  that  characterized  the  original  work.  The 
revised  model  retains  key  elements  of  the  original  model,  including:  ricochet  by  way  of  plastic  hinge  formation;  rod 
rebound  from  the  target’s  ricochet  surface;  and  target  gouging  as  a  means  to  redirect  the  target  force’s  line  of  action. 
Finally,  the  ricochet  equations  are  analytically  rendered  so  as  to  give  a  keen  sense  of  how  the  critical  obliquity  for  ricochet 
varies  as  a  function  of  the  initial  conditions  of  the  ballistic  engagement. 

©  2006  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

Recently,  Segletes  posited  a  model  for  rod  ricochet  [1,2]  based  on  the  notion  that,  during  the  ricochet 
event,  a  plastic  hinge  would  form  in  the  rod,  remaining  stationary  with  respect  to  the  target,  as  rod 
material  flowed  into  the  hinge  and  turned  through  an  angle  so  as  to  be  directed  away  from  a  penetrating 
aspect.  There  is  experimental  evidence  for  this  description  of  ricochet  morphology,  one  example  of  which  is 
shown  in  Fig.  1,  [3]. 

For  the  model  [1,2],  the  momentum  equations  for  rod  ricochet  were  derived,  based  upon  the  forces  and 
moments  acting  upon  and  fluxes  moving  through  a  control  volume  that  encompassed  the  plastic 
hinge,  as  described  in  Fig.  2.  By  assuming  that  the  incompressible  rod’s  cross-sectional  area  remains 
unchanged  through  the  course  of  the  ricochet,  continuity  dictates  that  the  rod’s  velocity  get  redirected 
but  remain  unchanged  in  magnitude  (F)  within  and  upon  exiting  the  plastic  hinge.  By  inferring  velocities 
directly  from  Fig.  1  measurements,  this  assumption  of  constant  V  can  be  shown  as  excellent  for  glancing 
obliquity.  Flowever,  its  validity  is  expected  to  diminish  as  the  striking  obliquity  becomes  more  normal. 
Nonetheless,  the  assumption  is  retained  for  all  impact  conditions  to  facilitate  analytical  solution.  Expres¬ 
sing  the  equations  for  linear  momentum  perpendicular  to  the  rod  axis  and  target-force’s  line  of  action. 
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Fig.  1.  Spark  cinematography  of  a  ricocheting  rod  projectile,  [3], 


respectively,  Segletes  obtained 

F  sin(0  +  t])  —  T  +  Vrh  cos(0  —  a),  (1) 

/  sin(0  -ft])  —  T  cos(0  +  t])+  Vm[cos{x  +  t?)  —  sin(0  +  t])],  (2) 

where  6  is  the  impact  obliquity  as  measured  from  the  target  normal,  a  is  the  angle  of  projectile  rebound  as 

measured  between  the  target  surface  and  the  trajectory  of  the  ricocheting  rod,  rj  is  the  angle  between  the  target 

normal  and  the  line  of  action  of  the  force  exerted  by  the  target  (positive  when  the  line  of  action  moves  away 
from  the  rod’s  initial  trajectory),  F  is  the  force  applied  by  the  target  to  the  plastic  hinge,  /  is  the  net  axial 
compression  force  in  the  rod  that  is  applied  to  the  hinge,  T  is  the  shear  traction  applied  by  the  cross-section  of 
the  rod  to  the  plastic  hinge,  and  Am/dt  is  the  rate  of  mass  flow  of  rod  material  through  the  plastic  hinge.  By 
assuming  a  frictionless  interface  between  rod  and  target,  17  =  0  for  those  cases  where  the  target  remains  elastic 
during  the  ricochet.  Only  when  the  target  is  plastically  gouged  during  ricochet  will  the  resultant  crater  allow 
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Fig.  2.  A  macro-view  of  ricochet  phenomenology,  depicting  forces  and  moments  upon  and  fluxes  through  the  plastic  hinge  contained 
within  the  shaded  control  volume  of  the  rod  [1], 


the  target  force’s  line  of  action  to  be  directed  at  a  nonzero  angle  t/  with  respect  to  the  target  normal.  In  essence, 
the  angle  rj  represents  an  angular  declination  of  the  effective  rod/target  bearing  surface  with  respect  to  the 
target  plane. 

The  notion  of  a  plastic  hinge  would  seem  to  limit  the  model’s  applicability  to  ductile  rods.  However,  the 
model  assumes,  as  revealed  in  Fig.  2,  that  that  ricocheting  material  exiting  the  hinge  region  provides  a  stress- 
free  boundary  condition  to  the  hinge’s  governing  momentum  equations.  In  this  sense,  regardless  of  whether 
the  ricocheting  material  is  intact  rod  (as  in  Fig.  1)  or  a  spray  of  redirected  fragments,  the  assumed  stress-free 
boundary  at  the  hinge’s  exit  seems  analytically  defensible.  And  while  the  mechanics  of  ductile  bending  and 
brittle  fracture  are  quite  different,  the  equations  being  solved  here  are  momentum  equations,  not  solid- 
mechanics  equations.  There  is  thus  reason  to  believe,  as  long  as  the  material  entering  the  hinge  can  support 
bending  stress  and  the  material  leaving  the  hinge  is  stress  free,  that  the  modeling  approach  may  retain  validity 
even  when  the  rod  is  brittle,  where  ricochet  is  accompanied  by  the  fracture  of  the  rod.  In  the  case  of  a 
ricocheting  spray  of  fragments,  the  model  variable  a  would  necessarily  refer  to  a  single  momentum-averaged 
trajectory  that  characterizes  the  spray. 

The  magnitude  of  the  momentum  flux  in  any  cross-section  of  the  rod,  is  given  as 

Vm^pV'^Ao,  (3) 

where  p  is  the  density  of  the  rod  material,  and  Aq  is  the  cross-sectional  area  of  the  rod.  The  maximum  axial 
force  that  the  rod  can  sustain  is 


/max  =  YA,,  (4) 

where  Y  is  the  strength  of  the  rod.  Segletes  introduced  an  “interaction  ellipse”  between  rod  and  target  by 
projecting  the  rod’s  cross-section  upon  the  target  surface.  By  further  introducing  the  artifice  of  a  flat  “bearing 
surface”  between  rod  and  target,  normal  to  the  action  of  the  target’s  applied  force  (i.e.,  canted  by  angle  f/  with 
respect  to  the  interaction  ellipse),  the  maximum  allowable  force  that  the  target  could  bear  upon  the  plastic 
hinge  of  the  rod  was  deduced  as 


F 


max 


^HAq 


cos  t] 
cos  9’ 


(5) 


where  FI  is  the  target’s  resistance  to  penetration. 
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Without  further  constraint,  there  is  the  possibility  of  many  solutions  to  the  momentum  equations,  Eqs.  (1) 
and  (2),  through  variation  of  the  rod  and  target  forces  over  the  allowable  domains.  Segletes  [1,2]  adopted  a 
moment-of-momentum  balance  as  an  additional  governing  equation  to  help  constrain  the  problem.  However, 
its  introduction  required  a  detailed  specification  of  the  stress  fields  acting  upon  the  plastic  hinge,  and 
complicated  the  solution  method  greatly.  Solutions  to  the  system  of  equations  were  obtained  through  a 
computational  approach  involving  significant  iteration.  In  an  effort  to  greatly  simplify  the  approach,  and 
remove  the  need  for  iterative  solution  methods,  we  herein  choose  to  dispense  with  the  moment-of-momentum 
balance,  and  replace  it  with  a  much  simpler  constraint  requiring  the  minimization  of  the  projectile’s  rebound 
angle  a. 

Such  a  minimization  can  be  justified  on  two  grounds:  (1)  the  rate  of  plastic  work  in  the  deforming  projectile 
is  minimized  (since  minimizing  the  rebound  angle  likewise  minimizes  the  plastic  turning  angle  of  the  rod);  and 
(2)  the  stress  load  upon  the  target  and,  by  inference,  the  target’s  stored  elastic  energy  are  minimized  during  the 
ricochet. 


2.  Revised  model  development 


Define  the  nondimensional  parameters 


A 


H 


pV^ 


(6) 


Y 

pV^ 


and 


(7) 


(8) 

The  parameters  A,  B,  and  /I  are  governed  by  the  material  properties  of  the  rod  (p,  Y),  the  target  (//), 
and  the  speed  of  the  impact  (V).  Note  that,  in  light  of  Eqs.  (3)-(5)  and  these  definitions, 
A  —  {FaiiixlVm)  ■  (cos  9/  cos  rj)  and  B  —  define  an  effective  value  of  A,  call  it  a, 

which  can  be  employed  to  represent  the  actual  target  force  during  ricochet,  F,  relative  to  the  maximum  allowed 
target  force,  as  a/ A  —  F/Fm^x-  Correspondingly  for  the  rod,  one  may  define  an  effective  value  of  B,  call  it  b,  in 
terms  of  the  actual  net  and  maximum  pennissible  rod  forces,  asb/B  —f  //max-  Likewise,  let  /Igjj-  —  \  b. 

For  cases  where  target  gouging  occurs  during  ricochet,  the  target  force  F  is  necessarily  at  its  maximum 
allowable  (i.e.,  plastic)  value,  iVax-  Thus, 

Target  force  constraint: 


0<a^A,  t]  —  0  (elastic  target), 
a  —  A,  t]^0  (plastic  target). 


The  rod,  however,  will  always  be  plastic  in  the  bending  hinge,  while  the  net  compressive  force  imparted  by 
the  rod  at  the  hinge  will  depend  on  the  location  of  the  neutral  stress  fiber  in  the  plastic  hinge.  For  example,  if 
the  neutral  fiber  of  the  hinge  exactly  traverses  the  centroid  of  the  rod’s  cross-section,  the  rod  will  be  in  a  state 
of  pure  bending,  and  /  will  be  zero  (even  as  the  rod  is  yielding,  bending,  and  ricocheting),  such  that  ^gff  will 
identically  equal  unity.  Constraining  the  rod  to  a  non-tensile  state  between  pure  bending  and  pure 
compression  gives: 

Rod-hinge  force  constraint: 

U/lgff^/?.  (10) 

The  original  momentum  equations,  Eqs.  (l)-(2),  when  shear  stresses  are  ignored  in  the  rod  (i.e.,  when 
T  =  0),  may  be  expressed,  via  direct  substitution,  as 
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Governing  momentum  equations: 

cos(9  —  a)  cos  9 

a  — - 

sin(0  +  q)  cos  q 

and 


(11) 


/leff  — 


cos(a  +  q) 
sin(0  +  q)  ■ 


(12) 


Failure  to  consider  the  possibility  of  shear  stresses  in  the  rod,  in  the  development  of  Eqs.  (11)-(12),  is  a 
regrettable  necessity  to  facilitate  the  solntion  pursned  in  this  article.  However,  a  mitigating  factor  in  this 
assnmption  is  that  the  shear  traction  T  must,  from  kinematic  arguments,  take  on  a  positive  magnitude,  using 
the  convention  of  Fig.  2.  A  positive  value  of  T  works  against  the  actions  of  both  Fin  Eq.  (1)  and /in  Eq.  (2), 
implying  that  target  stresses  and  the  axial  rod  force  must  increase  to  compensate  for  any  non-trivial  value  of  T. 
While  there  exists  one  very  limited  circumstance  in  which  this  sitnation  conld  be  helpfnl  (i.e.,  when  a  negative / 
wonld  otherwise  be  needed  to  facilitate  a  ricochet  solntion),  the  introduction  of  shear  traction  F,  in  general, 
works  against  the  viability  of  ricochet  solution. 

Qnite  simply  then,  with  the  omission  of  the  moment-of-momentum  equation,  the  remainder  of  this  report  is 
directed  to  the  analytical  solution  of  the  system  given  by  Eqs.  (11)-(12).  In  addition  to  the  force  constraints, 
Eqs.  (9)-(10),  a  list  of  kinematic  constraints  cnrrently  employed  is  given  below,  including  rebound 
minimization  and  those  developed  previously  [1,2]  that  still  apply  to  the  revised  approach. 


3.  Kinematic  ricochet  constraints 

(1)  Ricochet  definition  criterion:  ricochet  is  defined  as  when  the  trajectory  of  an  impacting  rod  rebounds 
away  from  (i.e.,  bounces  off)  the  surface  plane  of  the  target  or,  in  the  limit,  is  redirected  along  a  trajectory  of 
grazing  incidence.  Using  the  current  nomenclature,  a^O. 

(2)  Plastic-hinge  turning-limit  criterion:  a  tenet  of  the  ricochet  solution  is  that  the  rod  may  not,  by  way  of 
plastic  hinge,  turn  through  an  angle  greater  than  njl.  The  argnment  for  this  criterion  is  that  a  penetrator-dwell 
event  would  be  kinematically  more  likely  vis-a-vis  an  obtuse  flow-turning,  given  that  the  stresses  developed  in 
turning  a  rod  through  an  angle  greater  than  njl  wonld  exceed  the  stagnation  stresses  of  the  dwell  alternative. 
This  criterion  translates  to  the  statement  a^9.  Thus,  in  any  valid  ricochet  solution,  the  maximnm  permissible 
valne  of  a  is  0. 

(3)  Reflection-angle  limit  criterion:  given  that  q  represents  the  angular  declination  of  the  bearing  surface 
relative  to  the  target  plane,  the  criterion  that  limits  the  angle  of  reflection,  (a  -I-  q),  to  an  angle  less  than  or  eqnal 
the  angle  of  incidence,  {nl2—9—q),  follows  as  (x^nj2—9—2q.  This  criterion  snpersedes  the  one  posited  in  the 
prior  model  [1,2],  which  failed  to  account  for  the  declination  of  the  bearing  snrface  in  evalnating  the  angles  of 
incidence  and  reflection.  By  starting  with  Eq.  (12)  and  substituting  sin(7t/2— a— >;)  for  cos(a  -I-  q),  this  reflection- 
angle  limit  criterion  may  be  shown  eqnivalent  to  a  subset  of  Eq.  (10),  requiring  that  /ieff^l  (i.e.,  that  the 
impacting  rod  not  be  in  a  state  of  net  tension). 

(4)  Gouge-geometry  criterion:  based  on  arguments  [1,2]  nsing  Eig.  3,  it  was  shown  that,  for  any  symmetric 
ricochet  gonge  whose  lip  forms  an  angle  C  with  respect  to  the  target  surface,  the  line  of  action  of  the  target 
force,  defined  by  q,  and  the  rebound  angle  a,  were  constrained  by  the  gonge  geometry  in  such  a  way  that  \q\^oi 
is  a  requirement  (the  argument  conclnded  that  and  ^^\q\  simultaneously,  leading  to  the  constraint). 
While  solutions  with  negative  q  are  viable,  an  elastic  solntion  (in  which  q  =  0)  may  always  be  obtained  in  snch 
cases.  Thus,  for  practical  purposes,  q'^0  always  and  the  gonge-geometry  criterion,  in  application,  simplifies  to 
O^q^a. 

(5)  Rebound  minimization  criterion:  the  additional  constraint  which  will  serve  to  define  a  nniqne  solution  to 
the  momentnm  equations  (and  which  replaces  the  moment-of-momentum  constraint  ntilized  in  the  earlier 
model  [1,2])  is  that  the  rebonnd  angle  a  is  to  be  minimized  (subject  to  the  definition  of  ricochet  that  a^O).  As 
previously  stated,  this  minimization  constraint  has  the  effect  of  minimizing  both  the  rate  of  plastic  work  in  the 
ricocheting  rod  as  well  as  the  stored  elastic  energy  in  the  target,  dnring  the  event.  However,  and  like  the  prior 
methodology,  the  elastic  target  {q  =  0)  solntion  space  is  exhansted  before  considering  gouging  solntions 
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involving  non-trivial  values  of  r],  even  though  there  may  exist  certain  cases  in  which  a  gouging  solution 
produces  a  smaller  rebound  angle  than  the  corresponding  elastic-target  solution.  It  will  be  shown,  for  a  given 
choice  of  A  and  0,  that  to  choose  an  otherwise  viable  solution  with  the  smallest  a  is  equivalent  to  choosing  the 
viable  solution  with  the  largest 

4.  Elastic  ricochet 


The  case  of  elastic  ricochet  (i.e.,  for  which  the  target  remains  elastic)  presents  a  simplified  form  of  the 
general  ricochet  problem.  While  the  target  force,  embodied  in  variable  a,  is  known  in  advance  only  insofar  as 
the  bounds  0<a^^,  the  definitional  knowledge  that  r]  =  0  for  elastic  ricochet  greatly  simplifies  the  governing 
momentum  equations.  Derivations  for  the  elastic-  and  other  special-case  solutions  are  outlined  in  Appendix  A. 
Here,  only  final  results  are  presented.  These  results  follow  from  Eqs.  (9)-(12),  subject  to  the  various  kinematic 
constraints. 

The  basic  morphology  of  elastic  ricochet  is  that,  if  the  striking  obliquity  angle  0  is  large  enough,  the 
trajectory  of  the  incoming  rod  may  be  redirected  with  the  minimum  flow  turning  to  a  grazing  incidence 
parallel  with  the  surface  of  the  target  (i.e.,  there  is  no  rebound  of  the  ricocheted  rod  from  the  target  surface 
and  a  =  0  as  a  result).  In  these  cases,  the  rod  turns  through  an  angle  of  njl—O  within  the  plastic  hinge. 
Eventually,  as  the  striking  incidence  9  is  decreased  towards  the  normal,  there  comes  a  point  at  which  either  the 
target  strength  can  no  longer  sustain  the  increasing  load  associated  with  the  momentum  flux  of  the  incoming 
rod,  or  else  the  rod’s  strength  can  no  longer  sustain  the  large  momentum  flux  associated  with  the  ricocheted 
rod.  This  limiting  obliquity  at  which  zero  rebound,  a,  still  persists  is  designated  as  0c(  =  o- 

Segletes  [1,2],  in  an  early  attempt  to  simplify  the  calculation  of  ricochet,  actually  solved  and  presented  this 
very  restrictive  a  =  0  ricochet  result  as  a  “simplifled  model.”  Using  the  nomenclature  of  the  present  article,  the 
simplified  ricochet  criterion  was  given  as 


sin  0s,=o  =  max 


-A/2  + 


fl+iA/2f,l/p 


(13) 


For  some  cases,  however,  the  elastic-ricochet  event  can  be  prolonged  to  lower  9  if  the  penetrator  responds  by 
rebounding  off  the  target  surface.  The  rebound  angle,  a,  increases  as  9  is  successively  decreased  until  some 
critical  angle  of  elastic  ricochet,  6®,  below  which  elastic  ricochet  is  no  longer  possible.  At  this  critical  level  of 
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ricochet,  the  rebound  has  reached  its  maximum  elastic  value,  Fo'"  intermediate  values  of  impact  obliquity 
0l^6^6y  =  o,  the  rebound  falls  in  the  range  Altogether,  in  the  solution  of  the  equations  for  elastic 

ricochet,  there  are  three  branches  of  behavior  that  comprise  the  solution. 

In  what  is  called  “branch  1,”  it  is  the  target  strength,  embodied  in  A,  that  is  the  limiting  factor  which  prevents 
continued  ricochet  as  the  impact  obliquity  6  is  lowered.  Furthermore,  in  this  branch  of  ricochet  behavior,  the 
problem  cannot  be  remedied  by  way  of  an  increased  rebound  angle  a,  since  any  increase  in  a  comes  at  the  cost 
of  an  increased  target  stress.  For  branch  I  elastic  ricochet,  the  obliquity  limit  for  a  =  0  ricochet,  is  identical  to 
the  critical  angle  for  elastic  ricochet,  and  there  is  no  rebound  in  this  branch  of  elastic  ricochet: 

Elastic  ricochet  branch  I:  A<fi—  1  //i 


6^=0  =  sin“ 


-A/2  +  f\+{A/2f 


(14) 


01  =  0«=o,  (15) 

«Lx  =  0,  (16) 

At  0  =  0"  :  f  and  F  =  F^ax- 

At  the  other  extreme,  designated  “branch  III”  elastic  ricochet,  it  is  the  rod  strength  which  limits  the  non¬ 
rebounding  elastic  ricochet.  However,  in  this  solution  branch,  the  critical  ricochet  angle  may  be  lowered  below 
the  a  =  0  non-rebounding  limit,  thereby  inducing  the  rod  to  rebound  from  the  target  surface.  As  the  rod 
trajectory  changes  in  the  form  of  increasing  rebound,  in  response  to  a  lowering  of  the  striking  obliquity,  an 

increasing  fraction  of  the  rebounding  rod’s  momentum  flux  is  borne  as  load  by  the  elastic  target.  Assuming 

that  the  target  strength  can  elastically  sustain  this  increase  of  load,  this  compensating  process  can  continue 
with  decreasing  striking  obliquity  0  until  the  plastic-hinge  turning-limit  criterion  is  reached: 

Elastic  ricochet  branch  III:  A  >P 

0,=o  =  sin-‘(l/i3),  (17) 


0"  =  tan  '(1//1), 


(18) 


a 


e 

max 


^01 


(19) 


For  0^^0^0„=o  :  a  =  cos  HP  sin  0),  /  =/max  and  FF^ax- 

Between  these  two  extreme  branches  of  target-limited  and  rod-limited  ricochet,  “branch  11”  elastic  ricochet 
dehnes  the  region  where,  initially,  it  is  the  rod’s  strength  which  limits  the  obliquity  of  non-rebounding  ricochet. 
However,  as  the  striking  obliquity  is  further  lowered,  accompanied  by  increasing  rod  rebound,  the  increased 
target  load  associated  with  the  rebound  eventually  pushes  the  target  to  its  limit  force,  thereby  limiting  the 
obliquity  for  elastic  ricochet: 

Elastic  ricochet  branch  II:  P  —  l/ji^A^P 


0^=0  =  sin  il/P), 

(20) 

01  =  cos-‘  (yAj  f\+2A^-  , 

(21) 

“max  =  COS-'(^  sin  0^, 

(22) 

For  0l6^6a=o  :  a  =  cos”‘(^  sin  0),  f  and  FF^ax- 

At  0  =  0;:  :/  =/n,ax  and  F  =  F^ax- 

The  domain  of  engagement  scenarios  over  which  these  three  branches  of  elastic  ricochet  operate  is  depicted 
in  Fig.  4.  Therefore,  given  a  particular  ballistic  engagement  scenario,  one  may  ascertain  which  branch  of 
elastic  ricochet  is  potentially  operational.  Once  the  elastic-ricochet  branch  has  been  determined,  the  particular 
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Fig.  4.  Branches  of  elastic  ricochet,  depicted  as  a  function  of  rod  and  target  material  properties. 


values  for  00,  =  o  and  61  may  be  calculated  for  that  branch,  and  compared  with  the  actual  engagement  obliquity 
to  determine  whether  the  elastic-ricochet  conditions  are  met.  If  so,  all  that  is  left  is  to  determine  the  rebound 
angle  from  Eqs.  (16),  (19),  or  (22),  respectively.  The  calculation  is  complete  and  need  go  no  further.  Only  if  the 
engagement  obliquity  falls  below  the  critical  value  for  elastic  ricochet  need  further  consideration  be  given.  In 
this  case,  ricochet  involving  target  plasticity  (in  the  form  of  target  gouging)  needs  to  be  considered  as  a 
possibility.  The  next  section  on  plastic  ricochet  details  how  this  is  done. 


5.  Plastic  ricochet 


Below  are  given  the  solutions  to  various  special-case  conditions,  which  will  be  used  to  ascertain  the  more  general 
solution  for  plastic  ricochet.  For  all  special  cases,  one  may  assume  that  the  following  material  and  initial  conditions 
are  known  (i.e.,  fixed):  A,  B  (i.e.,  p  ),  and  6.  The  unknowns,  one  or  more  of  which  are  dehned  by  each  particular 
special  case  being  considered,  are  t],  a,  and  (S^ff  (i.e.,  b).  Because  the  target  is  known  to  be  plastic,  it  is  known  in 
advance  that  a  =  A.  These  special  case  solutions  are  not,  in  and  of  themselves,  guaranteed  to  satisfy  all  the 
necessary  constraints  of  ricochet.  However,  a  judicious  examination  of  the  constraints  in  concert  with  the  various 
special-case  solutions  will  permit,  if  it  exists,  a  general  solution  to  be  established  that  satishes  all  the  constraints.  To 
assist  the  reader  in  deriving  the  hnal  results  given  below,  derivation  strategies  for  these  special-case  solutions  are 
given  in  Appendix  A.  These  results  follow  from  Eqs.  (9)-(12),  subject  to  the  various  kinematic  constraints. 

5.7.  Solution  for  a  specified  fieff 


When  a  value  for  /ieff  is  directly  specified,  a  solution  for  the  governing  momentum  equations  may  be 
obtained  as 


'7/?eff  — 


1  sin"  ‘ 


where 


q^A^/2-  {^^f(A  +  1  -  ^eff)cos^  6, 


(23) 

(24) 


r  —  Afi  j 2  —  PgffA  cos^  6 


(25) 
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P  —  sill  ^  COS  9.  (26) 

Note  that,  while  negative  rj  solutions  are  valid,  they  will  never  prevail,  because  there  is  always  an  elastic  (rj  =  0) 
ricochet  solution  at  equal  or  lesser  value  of  /leff-  The  associated  value  of  a  is  given  as 

a/?erf  =  0-  sin”‘[jSeff  -  A  cos(0  +  cos  cos  0],  (27) 

5.2.  Solution  for  which  a.  =  0 

It  may  be  shown  that  an  a  =  0  solution  to  the  governing  momentum  equations  exists  if 


2  sin  9 


^  1  (if  a  =  0  solution  exists) 


is  satisfied  by  the  material  properties  and  initial  conditions  of  impact.  The  value  of  rj  when  Eq.  (28)  is  satisfied 
is  given  by 


'7«=o  =  2®^“^ 


1  .  _  I  2  cos^  9  sin  29 

sm  - ; - - — 

2  A  2 


1  +  1  / 1 


4  sin  9  4cosi  9 
'^4  A^ 


The  value  of  /leff  that  satisfies  this  condition  is 


Pa=o  — 


cos 

sin(0  +  ' 


Should  Eq.  (28)  remain  unsatisfied,  no  ricochet  is  possible,  as  it  can  be  shown  that,  for  such  a  condition, 
there  will  exist  no  solutions  for  which 

5.3.  Solution  that  maximizes  a 

Based  on  the  plastic-hinge  turning-limit  criterion,  the  maximum  permissible  value  for  a  is  9.  An  a  =  0 
solution,  which  represents  a^ax,  is  obtainable  if 


- - -  ^  1  (for  7.-9  solution  to  exist)  (31) 

4^  tan  0  ^ 

is  satisfied  by  the  initial  conditions  of  the  problem.  The  value  of  rj  when  Eq.  (31)  can  be  satisfied  is  given  by 


0a=0  =  2®^“^ 


f  .  _i  2cos^0  sin  20 


1  +  t/ 1  + 


4  tan  0  4 


The  value  of  /lefr  that  satisfies  this  condition  is 
^“=''  =  tan(0 +  ,?,=,)■ 

Should  Eq.  (31)  remain  unsatisfied,  however,  there  will  be  no  solution  for  which  a  =  9  .  Nonetheless,  a 
solution  may  be  obtained  which  maximizes  a  (while  keeping  «<  0)  by  minimizing  0— a.  In  such  cases,  it  may  be 
shown  that 

0(fl-c()min  =  ^/4  —  d/2.  (34) 

ft  is  useful  to  define  a  quantity  k  as  the  cosine  of  this  minimized  (0— a)  angle.  From  this  definition,  it  follows 
that 

Kmax  =  0  -  COS"‘  k,  (35) 

where  it  can  be  shown  that 

A(  l  +  sin  0) 
k  =  cos(0  -  a)„i„  =  ■ 

2  cos  0 


(36) 
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The  value  of  /Igff  satisfying  this  condition  is 
P(e-oi)mm  —  ~  +  ^/2- 


(37) 


5.4.  Solution  for  which  a  =  ri 

Based  on  the  gouge-geometry  criterion,  f/^a,  it  is  useful  to  know  the  conditions  that  provide  a  solution  to 
Eqs.  (11)-(12)  for  the  limiting  case  of  f;  =  a.  The  solution  for  this  case  may  be  reduced  to  a  fourth-order 
polynomial  in  tan  >7,  as 


\-A^ 


4(cos^  0  —  A^) 


f  1 


"  =  +  snre  •“  +  s.n  2fl  ^ 


(38) 


Standard  analytical  methods  [4]  may  be  applied  to  Eq.  (38),  so  as  to  obtain  the  roots  to  this  quartic 
polynomial  in  a  direct  fashion.  Of  the  four  solutions  obtained,  the  smallest,  positive,  real  solution  for  tan)7o,  =  ,, 
is  the  one  desired,  if  it  exists.  Thus, 

>7^^^  =  tan” '[smallest,  positive,  real  solution  to  Eq.  (38)]. 

While  solving  Eq.  (38)  for  the  case  of  very  large  A  values  (e.g.,  very  slow  impacts  which,  nonetheless,  deform 
the  target  plastically  while  forming  a  plastic  hinge  in  the  rod)  is  not  of  particular  ballistic  interest,  it  should  be 
noted  for  completeness  that  experience  reveals  that,  even  in  double-precision  arithmetic,  errors  of  precision  in 
the  analytical  solution  can  produce  unacceptable  results  in  which  valid  ricochet  scenarios  are  not  evaluated  as 
such.  It  may  be  tempting  simply  to  address  the  issue  by  realizing  that  valid  plastic-ricochet  scenarios  for  large 
A  can  only  occur  at  very  small  values  of  9  (large  9  situations  are  elastic)  and  that,  for  very  small  9,  the  value  of 
cos6~l.  In  such  cases,  from  Eq.  (II),  it  is  deduced  that  rj  is  nearly  independent  of  /?eff  (and  a),  with  a  limiting 
value  of 


^  sin-'(lM)-0.  (39) 

However,  a  slightly  more  accurate  approximate  solution  to  the  a  =  >7  problem  may  be  obtained,  solving  Eq. 
(11)  for  large  A,  by  first  assuming  that  a  may  be  approximated  by  t]*,  where  t]*  is  a  known  value  of  r]  is  the 
vicinity  of  the  a  =  ri  solution  being  sought.  Because  the  behavior  of  rj  is  nearly  constant  with  Petr  for  large  A, 
this  approximation  is  very  good.  Eliminating  a  thus  from  Eq.  (11)  permits  the  solution  to  be  reduced  from  a 
quartic  to  a  quadratic  equation,  the  result  being 


1 


-1 


tan  9 cos^  9(e  —  I  —  f  E{2  —  ii)tan^  9  +  1 


(for  large  A), 


where 


2  cos(0  —  >7*) 
A  tan  9 


(40) 


(41) 


In  the  actual  implementation  of  the  general  solution  algorithm,  when  this  large-^  difficulty  is  detected,  it  is 
convenient  to  use  the  value  rj*  —  tj^  though  this  is  by  no  means  the  only  possibility.  Once  ijaL  =  ,j  is  obtained 
thus,  a  is  reset  to  this  value  from  its  earlier  rj*  approximation. 

Einally,  regardless  of  the  method  by  which  >7a  =  ,j  is  obtained,  the  associated  fi  is  given  as 


cos(2tj^^^) 

■'>  ~  sin(0  -t-  ■ 


(42) 
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5.5.  Solution  for  which  t]  =  0 

In  a  like  manner,  the  momentum  equations  may  be  solved  for  the  special  condition  of  >7  =  0,  to  yield 

=  ^  +  V^l  -  M  tan  e)\  (43) 

a,,^o  —  0  —  cos“'(^  tan  6).  (44) 

Of  course,  it  follows  that  the  solution  for  this  condition  exists  only  when  A  tanO^  1. 

6.  General-solution  strategy 

The  special-case  plastic-ricochet  solutions,  Eqs.  (23)-(44),  show  that  there  are  many  solutions  which  satisfy 
Eqs.  (1 1)-(12).  And  so,  while  the  elastic-ricochet  solution  is  direct  and  straightforward  to  employ,  the  question 
naturally  arises,  in  connection  with  the  constraint  to  minimize  a,  how  to  select  the  proper  ricochet  solution 
when  the  target  is  no  longer  elastic. 

To  answer  this  question,  two  derivatives  need  to  be  established  (keeping  in  mind  that  A  and  6  are  fixed  for 
this  analysis).  Eirst,  take  Eq.  (11)  (which,  for  the  plastic  case,  a  =  A)  and  differentiate  with  respect  to  t]  to 
obtain 

^  ^  A  cosje  +  2ri) 
d>7  cos  9  sin(0  —  a)  ’ 

Over  the  “domain  of  viable  solutions,”  defined  here  as  the  solution  space  to  Eqs.  (1 1)-(12)  in  which  the  initial 
four  kinematic  constraints,  embodied  in  the  relation 

min  [0,  7i/2  —  0  —  2>7],  (46) 

hold,  it  is  clear  from  Eq.  (45),  in  light  of  Eq.  (46),  that  da/d>7>0.  This  means,  within  this  domain  of  viable 
plastic  solutions,  that  for  any  change  in  conditions  that  decreases  t],  a  will  likewise  decrease,  and  vice  versa. 
Next,  take  Eq.  (12)  and  differentiate  with  respect  to  Petr  to  obtain 


dry 


sin(0  -1-  rj) 


-  -  (47) 

d/^eff  /^eff  cos(0  +  t])  +  sin(a  -f  >y)(l  -t-  da/dry) 

Eikewise,  with  all  the  individual  terms  on  the  right  side  being  identifiable  as  positive,  within  the  viable  solution 
domain,  it  is  clear  that  d>y/d/ieff<0. 

What  can  be  drawn  from  Eqs.  (45)  and  (47)  is  that,  when  the  viable  solution  domain  (for  fixed  A  and  0)  is 
examined  as  a  function  of  /?eff,  both  rj  and  a  must  decrease  with  increasing  /?eff.  Therefore,  the  constraint  of 
minimizing  a  translates  into  one  of  selecting  the  maximum  value  of  flett  within  the  domain  of  viable  gouging- 
ricochet  solutions. 

For  the  case  where  an  a  =  0  solution  exists,  via  Eq.  (28),  for  a  given  A,  0  combination  (since  ricochet  is 
otherwise  precluded),  Fig.  5  describes  the  four  modes  of  behavior  for  a  and  >y  over  the  viable  solution  domain. 
The  figure  also  includes  a  very  rare  variant  of  Mode  IV  behavior  (deemed  Mode  IV-x)  in  which  a  double  a  =  >y 
crossover  occurs.  In  each  case,  different  logic  must  be  applied  in  order  to  obtain  the  largest  value  of  ^eff  in  the 
viable  solution  domain,  and  then  guarantee  that  the  solution  simultaneously  satisfies  the  rod-hinge  force 
constraint,  Eq.  (10). 

The  algorithmic  train  of  logic  that  will  accomplish  this  task  is  described  in  Table  1.  To  summarize  the  logic 
in  brief: 

If  no  elastic  solution  is  viable,  restrict  the  /ieff  domain  to  the  intersection  of  the  material  domain  [1,  /i]  with 

the  kinematic  domain  [/iamax.  Px  =  q]',  If  the  resultant  /igff  domain  is  not  the  null  set,  pick  the  maximum  /ieff  in 

this  domain  for  which  rj^a,  if  such  a  value  of  ySeff  exists. 

Note,  in  the  algorithm  (and  Fig.  5),  that  if  the  value  of  p,  associated  with  Kmax,  would  otherwise  fall 
in  the  invalid  [tensile]  domain  less  than  unity,  the  parameter  a^^^x  is  reset  to  the  value  associated  with 
^eff=l.  Furthermore,  in  the  actual  implementation  of  the  algorithm,  the  t]  =  0  special  solution  is  not 
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Mode  I;  7,  <  0, 


(a)  Arr 


Mode  II:  ^,  <  0,  rf^  > 


(b)  ^eff 


Mode  III:  >  0,  <  ^^,nax 


Mode  IV:  /7j  >  0,  //j  > 


Mode  IV-x:  >0,g^> 

(bul  with  double  a  =  g  crossover) 


(C) 


(d) 


(e) 


A„ 


Fig.  5.  Solution  to  the  plastic  ricochet  equations,  depicting  four  modes  (and  variant)  of  a,  rj  behavior  as  a  function  of  /?eff  (note:  iji  =  )7a_o)- 
Viable  solution  domain  arises  in  regions  where  the  kinematic  constraints  O^a^O,  /ieff>  1,  and  simultaneously  hold.  Because  it  is 

desired  to  minimize  a,  the  chosen  solution  will  be  the  largest  /?eff  in  the  viable  domain,  subject  to  the  rod-hinge  force  constraint,  /ieff  ^  P- 
Elastic  solutions  (shown  in  dotted  lines)  are  always  chosen  in  preference  to  plastic  ones. 


explicitly  used  to  limit  the  viable  solution  domain  for  /left,  because  it  is  implicitly  satisfied.  For  Modes  I 
and  II  behavior,  solutions  that  would  otherwise  satisfy  the  r]<0  condition  would  already  have  been 
discerned  as  elastic  solutions,  thus  precluding  the  plastic  algorithm’s  invocation.  For  Modes  III 
and  IV  behavior,  any  >7<0  solution  would  have  been  otherwise  ruled  out  by  the  limiter  associated 
with  the  a  =  0  criterion  which,  for  these  modes,  is  more  stringent  than  the  17  =  0  criterion.  Thus,  while  no  coding 
is  included  to  restrict  the  viable  solution  domain  to  positive  17  for  plastic  ricochet,  the  algorithm  nonetheless 
assures  it. 

Upon  ascertaining  the  values  of  a  and  r]  that  satisfy  the  particular  ricochet  solution,  the  force  quantities 
contributing  to  the  ricochet  are  also  of  interest.  Valid  for  situations  involving  either  elastic  or  plastic  ricochet, 
they  reduce  directly  from  the  momentum  equations,  given  the  definitions  of  a  and  b,  as 

F  a  \  cos(0  — a)cos0 

-fmax  A  A  sin(0  +  t])  cos  t] 

and 

/  hi/ cos(a  -f  >7)  ^ 


(49) 
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Table  1 

Ricochet  algorithm  logic 


Generate  elastic  ricochet  criteria,  as  a  function  of  given  A  and  fl. 

IF  0  greater  than  the  critical  angle  of  elastic  ricochet  (i.e.,  if  elastic  solution  viable),  THEN 
Obtain  elastic  ricochet  solution  for  ot  (keeping  r|  =  Q)^. 

ELSE 

Check  for  PLASTIC  ricochet  solution  as  follows: 

IF  A,  0  combination  does  not  permit  a  =  0  solutions,  THEN 

No  ricochet.  Without  a.  =  0  solution,  all  solutions  are  g<0  (not  valid). 

ELSE 

Initially  set  tentative  viable  solution  domain,  to 

Obtain  a  =  0  solution  (including  /i„-o  and  ))a_o).  Set  = 

Reduce  /imax,  if  necessary,  from  p  to 

Obtain  maximized  a  solution  [a  =  6  or  (0— a)n,ij. 

IF  (or  P^0-o:)mm)  larger  than  THEN 

Raise  from  1  to  /?x_o  (or  P(o-a)min),  accordingly. 

Set  1]  2  to  value  of  (or  i7(o_x)min)- 
ELSE 

Obtain  solution  for  specified  /ieff  equal  to  p^j^  =  1 . 

Set  ij2  to  value  of  Reset  o;,nax  to  value  of  ayjmin. 

END  IF 

IF  there  still  exists  a  nontrivial  domain  of  viable  p^rr,  [Pmin,Pmax],  THEN 
Obtain  largest  p^{[  in  domain  lPminfmsx]foi'  which  tts  follows: 

IE  (ji  ^  0  (Modes  I  or  II),  THEN 

Obtain  tentative  ricochet  solution,  by  specifying  P^jf  as  /imax- 
IF  ijeffXZeff  (i  e.,  P  outside  viable  Mode  II  domain),  THEN 

No  ricochet.  Max  allowed  p  on  wrong  side  of  ix,rj  crossover. 

END  IF 

ELSE  IF  )/2<s(max  (Mode  III),  THEN 

Obtain  tentative  ricochet  solution  at  a  = crossover. 

IF  (i.e.,  p  in  middle  of  viable  domain),  THEN 

Specify  PeS[=P  ricochet  solution,  to  restrain  Pew^p. 

END  IF 
ELSE 

IF  aL  =  r\  crossover  nonetheless  exists,  THEN 
l¥  Px-,,^P,  THEN 

Mode  IV-x  solution  at  /ieff=fe-n- 

ELSE 

IE  rtf,fia.p,  THEN 

Mode  IV-x  solution  at  p^[f=p. 

ELSE 

No  ricochet.  Mode  IV-x  value  of  P  too  low. 

END  IF 

END  IF 

END  IF 

No  ricochet.  Mode  IV:  g  everywhere  greater  than  a  over  domain  of  viable  p. 
END  IF 
ELSE 

No  ricochet.  Effective  domain  of  viable  P  nonexistent. 

END  IF 
END  IF 
END  IF 


“Critical  ricochet-decision  points  are  given  in  underlined  text. 


While  it  is  not  shown  here,  the  same  technique  employed  by  Segletes  [1,2]  may  be  used  to  relate  the  relative 
rod  force  ///max  to  the  location  of  the  neutral  fiber  in  the  plastic  hinge  and  hence  to  the  strain  field  in  the  rod. 
Such  knowledge,  along  with  the  rod’s  material  characteristics,  may  inform  an  estimate  of  whether  the 
ricocheting  rod  retains  its  integrity  or  breaks  into  a  spray  of  ricocheting  fragments. 
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1.  Comparison  with  prior  model  formulation 

It  has  been  noted  that  the  difference  between  the  model  derived  herein  and  its  predecessor  [1,2]  is  in  the 
replacement  of  a  moment-of-momentum  relation  with  a  rebound-minimization  constraint,  adopted  here  so  as 
to  simplify  the  model  and  afford  analytical  solution  to  the  governing  equations. 

Despite  what  might  seem  a  significant  change  to  the  model,  the  reality  is  that  the  identical  momentum 
equations,  solved  in  both  models,  serve  as  the  primary  determining  factor  of  ricochet  behavior.  One  might 
expect,  since  the  moment-of-momentum  relation  is  a  more  limiting  “equality  constraint”  than  the 
“minimization  constraint”  that  replaced  it,  that  the  current  model  permits  a  slightly  larger  solution  space 
for  ricochet.  We  see  in  Fig.  6,  which  depicts  a  model  comparison  by  way  of  a  ricochet  phase  diagram,  that  this 
is  true.  The  engagement  scenario  depicted  in  this  figure  was  a  test  case  presented  in  the  original  report,  in 
which  a  rod  of  density  2700  kg/m^  and  strength  T=1.2GPa  was  impacted  against  a  target  of  resistance 
H=  1.25  GPa.  Since  the  propensity  to  ricochet,  as  measured  experimentally  for  other  scenarios  and  reported  in 
the  prior  work,  was  at  least  as  great  as  any  of  the  ricochet  models  (and  frequently  greater),  the  expanded 
ricochet  domain  of  the  current  model  is  a  further  step  in  the  right  direction. 

Additionally,  since  the  current  model  employs  a  constraint  that  expressly  adopts  the  solution  which 
minimizes  projectile  rebound,  one  might  anticipate  that  the  current  model  will  tend  to  produce  ricochets  with 
smaller  values  of  a  than  the  original  model  [1,2].  This  too  has  been  verified.  For  the  simulation  of  Fig.  6,  when 
the  striking  velocity  is  selected  as  F=  700  m/s,  the  prior  work  reported  modeled  rebound  angles  approaching 
28°,  with  the  corresponding  fy  values  between  0°  and  11°.  In  contrast,  the  current  model  for  the  same 
simulation  parameters  places  both  a  and  tj  between  0°  and  14°  over  the  domain  of  ricochet  solutions. 

For  the  two  cases  that  were  compared  against  experimental  results  in  the  original  reports  [1,2],  the  revised 
formulation  of  the  current  model  predicts  ricochet  characteristics  identical  to  the  simple  model  of  ARL-TR- 
3257  [1],  for  the  case  of  the  short  LjD,  brittle  penetrator  (B32),  over  the  range  of  experimental  data  (Fig.  7). 
For  the  longer,  L/Z)=15,  93%  tungsten  rod  results,  critical  ricochet  obliquities  from  the  currently  revised 
formulation  are  2-6°  lower  than  even  the  simplihed  form  of  the  earlier  model  over  the  striking  velocity  range 
of  600-1200  m/s  (where  ricochet  was  experimentally  observed),  bnt  still  not  quite  to  the  critical  obliquity  level 
noted  in  the  experiments  (Fig.  8).  The  original  work  [1,2]  may  be  consulted  for  specific  details  of  the  geometry 
and  material  characteristics  of  the  engagements.  In  any  event,  the  revised  model  presented  here  produces  the 
most  accurate  prediction  of  ricochet  of  the  several  models  studied,  including  those  by  Tate  [5]  and  Rosenberg 
et  al.  [6]. 


Fig.  6.  Comparison  of  ricochet  phase  diagrams  for  engagement  characterized  by:  p  =  2700kg/m^,  Y=  1.2  GPa,  H=  1.25  GPa. 
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Fig.  7.  Comparison  of  ricochet  models  for  B32  penetration  into  aluminum  (note:  Tate  and  Rosenberg  models  predict  no  ricochet). 


Fig.  8.  Comparison  of  ricochet  models  to  L/T)=  15  tungsten  alloy  rods  into  aluminum. 


8.  Results  of  interest 

In  addition  to  comparing  the  revised  formulation  to  the  earlier  model,  it  is  interesting  to  explore  the  new 
formulation  for  its  own  sake.  One  point  that  should  be  noted  is  that,  while  most  of  the  equations  have  been 
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Fig.  9.  Sample  ricochet  phase  diagram  for  revised  model. 


devoted  to  solving  the  general  problem  of  gouging  ricochet,  only  a  small  portion  of  the  V,9  engagement 
envelope  is  characterized  by  it.  Fig.  9  shows  this  point  clearly  for  a  typical  steel-on-steel  engagement,  where 
much  of  the  diagram  is  either  elastic  a  =  0  ricochet  or,  alternately,  no-ricochet. 

Also  added  on  the  graph  is  the  reevalnated  ricochet  dividing  line  for  a  hypothetical  case  of  a  much 
stronger  rod,  which  is  shown  to  make  several  instrnctive  points.  First,  the  kink  in  the  F=lGPa 
result  at  around  350  m/s  impact  velocity  is  clearly  related  to  a  limitation  of  the  rod’s  strength  to  snstain 
ricochet,  as  this  kink  is  absent  for  the  “large  Y  case”  curve.  In  fact,  when  ricochet  occurs  at  low  striking 
obliquities,  the  rod  must  necessarily  turn  approximately  90°  during  the  ricochet,  and  travel  with  a 
correspondingly  small  rebound  from  the  target  snrface.  This  means  that  the  portion  of  rod  momentum  coming 
in  to  the  plastic  hinge  can  contribute  virtually  no  force  to  balance  the  momentnm  flux  leaving  the  hinge. 
Likewise,  the  target’s  force  is  nearly  perpendicular  to  the  flow  of  the  ricocheted  rod  and  can  also  contribnte 
only  a  small  component  of  its  force  to  balance  the  ricocheted  momentum  flux.  Therefore,  a  primary 
contributor  to  that  flux  balance  (and  what  therefore  is  a  limiting  factor)  must  necessarily  be  the  strength  of  the 
rod  applied  to  the  hinge. 

While  this  explanation  is  supported  by  the  equations,  the  predictions  of  ricochet  at  very  low  striking 
obliquities  must  be  taken  with  some  reservations,  since  they  predict  that  a  strengthened  rod  will  be  more  likely 
to  ricochet.  Consider  first,  that  the  alternative  to  ricochet  is  not  solely  penetration,  but  might  also  include 
dwell  (when  the  target  resistance  H  exceeds  the  dynamic  stress  of  the  impacting  rod,  Y+ljlpV^).  In  the 
example  of  Fig.  9,  dwell  would  occur  to  a  point  beyond  500  m/s  impact  speed  and  so  hypothetically 
strengthening  the  rod  does  not  change  the  scenario  from  penetration  to  ricochet,  but  from  dwell  to  ricochet, 
which  is  much  more  believable.  Also,  the  whole  basis  of  the  present  model  is  one  that  evalnates  whether  the 
rod  and  target  strengths  can  sustain  ricochet,  not  whether  they  can  initiate  it.  Becanse  it  is  known  that  ricochet 
is  very  nose-shape  dependent,  and  that  ricochet  can  not  be  snstained  if  it  fails  to  initiate  in  the  first  place,  the 
current  model’s  nnderlying  premise  makes  it  ill  snited  to  predict  whether  a  strengthened  rod  wonld  have 
initiated  ricochet  in  the  first  place. 

Another  interesting  behavior  was  observed  that,  at  first  glance,  might  seem  counterintuitive.  Recall  from 
Fig.  5  and  the  associated  discnssion  that  a  and  rj  were  determined  to  vary  monotonically  with  increasing 
Snch  monotonicity  in  /?eff,  however,  implies  no  such  comparable  behavior  as  a  function  of  6.  Using  the  same 
material  parameters  of  Fig.  9  and  considering  a  specific  impact  speed  of  F=  400  m/s,  the  variation  of  a  and  rj 
are  presented  as  a  function  of  9  in  Fig.  10.  Recall  that  these  curves  do  not  represent  pnrported  behavior  dnring 
a  single  impact,  but  that  each  chosen  obliqnity  9  represents  a  unique  ballistic  engagement. 
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Fig.  10.  Ricochet  response  of  a  and  jj  as  a  function  of  striking  obliquity. 


While  the  fluctuation  of  a  may  seem  erratic,  each  segment  of  behavior  is  understandable  in  terms  of  the 
constraints.  As  the  striking  obliquity  is  lowered  to  33.9°,  the  rod  strength  becomes  unable  to  sustain  a  =  0 
ricochet  and  starts  to  compensate  with  rod  rebound  off  the  elastic  target.  However,  not  only  does  lowering  0 
add  load  directly  upon  the  target  from  the  impacting  portion  of  the  rod,  but  the  increasing  rebound  of  the 
ricocheting  rod  further  increases  the  target’s  load.  When  the  target  strength  can  longer  sustain  this  load 
elastically  (at  a  striking  obliquity  of  31.6°),  it  begins  to  gouge  at  impact. 

The  gouging  allows  the  target  to  redirect  its  line  of  action  by  a  few  degrees  in  the  direction  of  the  ricocheting 
rod.  This  redirection  of  target  force  (an  increasing  rj)  provides  a  force  component  parallel  with  the  target 
surface  that  allows  the  momentum  flux  of  the  ricocheting  rod  to  be  sustained  at  a  lower  value  of  rebound  angle 
o£.  This  mode  of  compensation  (increasing  rj  in  exchange  for  decreasing  a)  proceeds  as  the  striking  obliquity  is 
lowered,  until  the  point  at  25.7°  striking  obliquity  where  the  value  of  r]  rises  to  meet  that  of  the  lowering  a.  At 
this  point,  the  gouge  lip  has  formed  a  7°  angle  ^  with  respect  to  the  target  surface  (see  Fig.  3).  The  target-force 
redirection  rj  must  be  at  an  angle  no  greater  than  this  gouge  angle,  while  the  rod  rebound  angle  a  must 
kinematically  be  at  least  as  large  as  Thus,  as  the  striking  obliquity  is  further  lowered,  the  target  gouge 
deepens,  which  allows  for  an  increasing  t]  but  only  at  the  kinematic  cost  of  increasing  the  rebound  angle  such 
that  is  maintained. 

Finally,  as  the  obliquity  is  lowered  to  19.45°,  the  rebound  a,  being  forced  upward  by  the  deepening  gouge, 
reaches  this  same  level  of  19.45°.  At  this  point,  the  rod  is  making  a  90°  turn  in  direction  from  its  striking 
orientation.  The  only  way  in  which  ricochet  could  be  sustained  at  lower  striking  obliquity  would  be  for  the  rod 
to  continue  turning  in  excess  of  90°.  However,  such  an  adjustment  is  considered  untenable,  since  the  forces 
required  to  achieve  this  would  exceed  the  forces  required  to  dwell  the  penetrator.  Thus,  for  impact  obliquities 
below  19.45°,  the  more  plausible  scenario  is  that  the  rod  would  simply  flatten  and  dwell  upon  the  surface  of  the 
target,  in  the  mode  of  a  Taylor  impact  cylinder. 

9.  Critical  ricochet  obliquity  relationships 

To  this  point,  all  calculations,  including  those  that  portrayed  a  critical  obliquity  curve  in  one  form  or 
another,  as  in  Figs.  6-9,  were  done  by  systematically  applying  the  aforementioned  general-solution  strategy 
over  the  range  of  initial  conditions  considered.  Such  a  method  is  necessary,  if  detailed  ricochet  response,  such 
as  that  portrayed  in  Fig.  10,  is  sought.  However,  when  all  that  is  required  is  the  minimum  obliquity  at  which 
ricochet  can  possibly  occur  (i.e.,  the  critical  ricochet  obliquity)  for  a  specified  combination  of  rod  and  target 
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material  descriptions  and  impact  velocity,  it  is  desirable  to  have  simpler  relationships  for  making  a  quick 
assessment  of  ricochet  likelihood.  The  development  of  such  relationships  would  also  serve  to  provide  a 
heightened  understanding  of  the  ricochet  phenomenon. 

When  the  rod  strength  is  large,  it  has  been  observed  in  the  model  that  the  critical  angle  of  ricochet  for  low 
striking  obliquity  occurs  when  6  =  x  =  t].  Since  the  kinematic  constraints  require  that  and  a^6,  this 
9  =  a  =  i'i  condition  represents  a  clear  limiting  case  for  the  critical  ricochet  angle.  Based  on  the  kinematic 
constraint  that  however,  this  condition  can  only  occur  at  striking  obliquities  below  9^22.5°. 

Thus, 

=  9c  (when  p  ^  V -  1,  0  ^  0c  ^  22.5°),  (50) 

where  the  /?  specification  quantifies  the  rod-strength  requirement.  The  critical  obliquity  for  these  conditions, 
from  the  momentum  equations,  may  be  derived  as 

9c  —  ^sin“'(l/^)  (when  P^Va^  —  1,  A^V2)  (51) 

or,  conversely, 

(when  -  1,  0^0^22.5°).  (52) 

sin  29c 

It  is  the  simultaneous  solution  of  Eq.  (12),  in  light  of  Eq.  (50)  (i.e.,  l/tan20c),  with  Eq.  (52)  that  provides 
the  constraint  between  A  and  /?.  The  \/2  limit  on  A  corresponds  to  that  value  which  will  bring  9c  to  22.5°.  It  is 
exactly  this  curve  which  corresponds  to  the  large- T  curve  extension  presented  in  Fig.  9.  Correspondingly,  it 
has  been  noted  in  the  model  that  if  the  /?  constraint  can  not  be  satisfied  for  a  given  A  in  the  range  A  >  ~/2,  then 
there  are  no  plastic  ricochet  solutions,  since  all  would-be  plastic  solutions  violate  the  plastic-hinge  turning- 
limit  criterion.  There  will,  however,  be  branch  III  (i.e.,  A>fi)  elaslic  solutions  for  A>V2,  should  the  /I 
requirement  not  be  met  for  providing  a  plastic  solution.  However,  these  elastic  solutions  occur  at  higher  values 
of  striking  obliquity,  resulting  in  a  jump  in  the  9c  curve. 

In  an  adjacent  range  of  critical  plastic-ricochet  obliquities,  immediately  above  22.5°  for  which  A<V2,  the 
limiting  ricochet  is  governed  by  the  gouge-geometry  constraint  and  the  reflection-angle-limit  constraint  which, 
taken  at  their  limiting  values,  give  Xc  =  r]c  and  ac  =  7i/2— 0c— 2f/c-  The  latter  relation,  in  light  of  the  former, 
reduces  to 

otc^rjc^  30°  -  0c/3.  (53) 

Furthermore,  the  reflection-angle  limit  constraint  indicates  that  this  condition  must  occur  at  a  value  of 
/leff=  1.  These  values  of  rjc  and  0c  are  then  substituted  into  Eq.  (11)  in  order  to  get  the  A  value  associated 
with  that  critical  obliquity: 

^  _  cos((40c/3)  -  30°)  cos  0c 

“  sin((20c/3)-t-3O°)cos((0c/3)- 30°)  ■  ' 

But  through  what  obliquity  beyond  22.5°  will  this  condition  prevail?  From  the  aforementioned  general- 
solution  strategy,  it  is  known  that,  when  the  solution  departs  from  this  condition,  /leff  will  increase  beyond  a 
value  of  unity,  and  Kc  will  diminish  in  comparison  to  the  [/leff=  1]  value  given  by  Eq.  (53).  Therefore,  we 
conduct  a  perturbation  analysis,  in  which  we  hypothesize  a  departure  from  this  solution  by  allowing  otc  to 
diminish  by  an  infinitesimal  amount  e: 


Kc  =  'Ic  =  30°  -  0c/3  -  £, 

(55) 

leading  to 

^  cos((40c/3)  —  30°  -t-  £)  cos  0c 

“  sin((20c/3)  +  30°  -  s)  cos((0c/3  -  30°  -t-  s)  ' 

(56) 

In  order  for  this  perturbation  from  the  baseline  solution  to  feasibly  establish  a  lower  critical  ricochet  angle 
(and  thus  an  alternative  to  the  baseline),  00c/0e<O  is  required.  Expressing  this  as  (0H/0£)/(0^/00c)<O,  we  note 
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that  0^/06c  will  be  negative  for  all  0^  in  the  range  22.5°  <  9c^90°,  and  so  the  feasibility  of  a  departure  from  the 


1  solution  hinges  on  the  sign  of  dA/dE  as  £->0.  It  may  be  shown  that  0^/0e  changes  sign  in  the  limit 
(£->0)  when  the  condition 


(57) 


is  met.  This  condition  is  met  when  0^  =  38.06°.  Thus,  the  range  on  this  ^eff=  1  solution  spans 
22.5°  <  0c ^38.06°,  through  which  ac  =  f;c  =  30°— 0c/3.  From  Eq.  (54),  it  may  be  observed  that  the  range  on 
A  for  this  region  is  0.9373 \/2. 

At  larger  0c,  all  the  way  to  90°,  there  continues  to  be  a  solution  governed  by  the  constraint  (Xc  =  ?7c  =  30°— 0c/ 
3,  for  /ieff=l,  though  this  particular  solution  is  no  longer  that  which  minimizes  a  for  cases  of  larger  /i. 
Nonetheless,  for  this  very  special  /ieff=  1  case,  Eq.  (54)  applies  for  relating  A  to  0c.  Further,  it  allows 
generalization  over  the  range  38.06° < 0c ^42.497°,  corresponding  to  the  range  0.8047^ A <0.9373.  In  this 
narrow  range,  for  the  A,  6^  pair  given  by  the  Eq.  (54),  there  are  two  positive-?^  solutions  for  the  a  =  r] 
constrained  momentum  equations,  i.e.,  for  Eq.  (38):  the  one  given  at  /ieff=  1,  and  the  other  with  ^gff  in  the 
range  between  1  and  1 .4803 .  For  this  second  solution,  the  resulting  value  of  =  9c  is  smaller  than  the  value  for 
the  /ieff=  1  solution,  and  is  thus  the  preferable  solution  for  ricochet. 

It  would  be  convenient  to  hope,  for  a  given  value  of  A,  that  when  0^  is  lowered  to  the  point  where  the  /?  =  1 
solution  disappears,  that  all  the  /ieff>  1  solutions  would  simultaneously  disappear.  Were  that  the  case,  then  the 
/i  =  1  solution,  while  not  the  preferable  solution  for  ricochet,  would  still  be  indicative  of  the  critical  ricochet 
angle,  0c.  However,  there  is  a  slight  imperfection  with  such  a  hypothesis,  and  that  is  the  phenomenon  of  Mode 
IV-x  ricochet  depicted  in  Fig.  5e,  in  which  the  P=  \  solution  is  no  longer  active,  despite  a  range  of  intermediate 
/i  for  which  a  ricochet  solution  does  exist.  This  mode  of  solution  can  exist  for  striking  obliquities  greater  than 
38.06°,  and  will  be  discussed  in  more  detail  later.  Suffice  it  to  say,  at  this  point,  that  the  region  where  Mode  IV- 
X  ricochet  can  operate  is  a  miniscule  triangular  region  of  the  A, 0c  plane.  Therefore,  while  not  exact,  Eq.  (54), 
relating  A  to  0c  for  the  ^  =  1  solution,  provides  an  excellent  approximation  to  the  critical  obliquity  for 
0.8047^ A <0.9373,  for  all  ^3^1. 

A  further  difference  between  this  solution  and  that  for  22.5° <  0c  =^38.06°  is  in  the  relationship  for  f/c  and  o£c. 
When  P=\,  the  solution  given  by  ac  =  >7c  =  30°— 0c/3  still  applies.  However,  when  /i  is  large  enough,  the 
preferable  secondary  solution  applies,  producing  a  lower  magnitude  of  ac  =  0c-  When  the  value  of  /i  is  in  a 
range  between  these  two  solutions,  the  value  for  will  be  less  than  that  for  Kc,  yet  both  will  fall  in  a  range 
between  the  two  oic  =  9c  solutions.  In  either  case,  is  diminished  with  respect  to  the  value  associated  with  the 
/i  =  1  solution. 

The  upper  end  of  this  range,  associated  with  the  A,  9^  pair  of  (0.8047,  42.497°),  corresponds  to  that 
condition  for  which  the  secondary  =  rjc  solution  is  exactly  ac  =  f^c  =  0,  traditionally  associated  with  an  elastic 
Branch  I  solution.  In  fact,  the  upper  end  of  the  range  was  ascertained  by  equating  the  value  of  A  for  the  /ieff  =  1 
plastic  solution,  Eq.  (54),  to  that  from  the  limiting  (a^A)  elastic  Branch  I  solution,  given  by  Eq.  (11)  with 
a  =  t]  =  0  and  a^A,  such  that 


cos^  0c  cos((40c/3)  —  30°)  cos  0c 


sin  0c  sin((20c/3)-t-3O°)cos((0c/3)-3O°)'  ^ 

It  is  the  solution  of  this  equation  which  yields  0c  =  42.497°,  corresponding  to  A  =  0.8047.  Thus,  it  is  not 
surprising  to  suspect  that,  at  obliquities  in  excess  of  42.497°,  elastic  solutions  may  influence  the  critical 
obliquity. 

Indeed,  when  ignoring  the  Mode  IV-x  solutions  in  the  range  of  critical  obliquities  0c  ^42.497°,  the  elastic 
Branch  I  criterion  may  govern  the  ricochet.  When  fi  exceeds  the  limit  needed  to  force  branch  I  elastic  ricochet 
(a  =  0)  in  this  range,  by  satisfying  the  relation  0.8047 </i—l//i  to  yield  /i>  1.4803,  Eq.  (11)  becomes 


C 


sin  0' 


(59) 
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Table  2 

Relations  for  critical  ricochet  obliquity 

A,  Oe  Range 

Relation^ 

Restriction 

ly,  a  Behavior  at  0c 

0c^22.5°  plastic, 

^=-\o 

Sin  20c 

P>Va^  -  1 

=  Oc  = 

0o^45°  elastic 

Elastic:  /?  =  — 

tan  01 

l^/tscVri^  -  1 

1]  =  0,  (x  =  6c 

0.9373^.4<2‘'^ 

22.5°<0oS:38.06° 

^  cos((40c/3)  —  30°)  cos  0c 

“  sin((20c/3)  +  30°)  cos((0c73)  -  30°) 

all  /?>! 

ac  =  ;7c  =  3O°-0c/3 

0.8047^yl<0.9373 

38.06°  <00^42.497° 

^  cos((40c/3)  —  30°)  cos  0,. 

~  sin((20c/3)  -1-  30°)  cos((0c/3)  -  30°) 

all  /?>! 

0^r]c^iXc^30'^—9J3 

0<ri<0.8047 

0e>42.497° 

^  cos((40c/3)  —  30°)  cos  0c 

'  “  sin((20c/3)  -1-  30°)  cos((0c/3)  -  30°) 

IS=l 

ac  =  J7c  =  3O°-0c/3 

.42^ri^.4i 

l^/lsS  1.4803 

0  ^  17c  ^  Kc  ^  30°— 0c/3 

cos^  0^ 

Elastic:  A2  ^  — itr  (relation  is  exact  for  0c  >45°) 

sin  0c 

1.4803 

ac  =  i7o~0 

^Relation  specified  is  for  plastic  ricochet,  unless  otherwise  noted.  Also,  the  presence  of  Mode  IV-x  ricochet  solutions  makes  some  of  the 
specified  relations  only  approximate,  denoted  by  the  symbol  a; ,  instead  of  = . 


Eq.  (59),  of  course,  can  be  inverted  to  yield  the  associated  elastic-ricochet  solution  of  Eq.  (14): 

0«  =  sin-‘  -A/1+  (when  iS>  1.4803,  0c  >42.497°  [i.e.,  when  0^y4<  0.8047]).  (60) 


Once  again  ignoring  the  influence  of  Mode  IV-x  solutions,  all  solutions  in  this  range  will  be  bounded,  as  a 
function  of  between  this  elastic  /!>  1.4803  solution  and  the  /?=  1  plastic  solution  of  Eqs.  (53)-(54).  Likewise, 
for  this  range  of  0^,  the  values  of  and  rj^.  will  range  between  the  values  associated  with  the  elastic  and  the 
/I  =  1  plastic  solution.  These  elastic  and  plastic  critical-obliquity  functions  are  actually  quite  narrowly  banded, 
producing  critical  obliquities  within  3°  of  each  other,  for  all  ^<0.8047. 

All  of  these  relations  described  for  critical  ricochet  are  collated  in  Table  2,  for  convenience. 

We  return  briefly  to  the  situation  of  the  Mode  IV-x  ricochet  solutions  which,  when  ignored, 
makes  the  described  critical  obliquity  description  an  inexact  one.  The  question  to  be  answered  is 
one  of  “how  inexact?”  Because  the  solution  to  the  a  =  r]  condition,  on  which  the  Mode  IV-x  solutions 
depend,  is  a  fourth-order  equation,  no  simple  description  has  yet  been  developed  for  the  precise  range  over 
which  this  mode  of  solution  may  apply.  However,  through  a  systematic  approach,  the  range  of  Mode  IV-x 
solutions  were  obtained,  in  order  to  show  the  small  degree  to  which  they  depart  from  the  analytical 
descriptions  otherwise  put  forth. 

Fig.  1 1  documents  the  result.  In  the  figure,  solid  lines  of  the  graph  denote  the  analytical  functions  presented 
here,  while  all  dotted  and/or  dashed  lines  represent  calculated  quantities  obtained  through  a  systematic 
solution  to  the  ricochet  equations  for  specifled  values  of  6,  A,  and  /I.  It  is  easy  to  see  how  good  a  job  the 
analytical  description  of  critical  ricochet  obliquity  angle  does,  even  allowing  for  the  Mode  IV-x  solutions.  In 
all  cases,  the  error  in  the  predicted  value  of  A  for  a  given  6^  is  always  under  1.2%,  and  generally  well  under 
1%.  Mode  IV-x  ricochet  plays  a  role  in  the  critical  obliquity  only  within  the  small  triangular  (A, 6 a)  region 
delimited  by  the  vertices  (1/V2,  45°),  (0.9373,  38.06°),  and  (0.8047,  42.497°).  Further  discussion  on  the 
intricacies  of  Mode  IV-x  ricochet  is  offered  in  Appendix  B. 

Fig.  12  summarizes  the  discussion  of  critical  ricochet  obliquity  as  a  function  of  target  parameter  A.  Note 
that  the  inclusion  of  plastic  and  elastic  ricochet  curves  in  this  flgure  is  not  meant  to  denote  a  transition  from 
one  ricochet  mode  to  another  with  increasing  obliquity  (though  such  transitions  indeed  happen).  Rather,  the 
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Fig.  11.  The  influence  of  Mode  IV-x  solutions  on  the  calculation  of  critical  ricochet  obliquity.  Solid  lines  represent  the  analytical  functions 
given  for  plastic  and  elastic  ricochet.  Actual  general-model  behavior  is  depicted  with  dashed  lines. 
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1/A'“  =  V/(Hlpf‘ 


Fig.  12.  Critical  ricochet  obliquity  as  a  function  of  parameter  A,  showing  the  interplay  of  plastic  and  elastic  solutions,  based  on  value  of  fl. 
Results  are  shown  (a)  for  small  values  of  A,  and  (b)  as  function  of  1/A'^^. 


multiplicity  of  curves  is  meant  to  show  how  the  critical  obliquity  for  ricochet  can  change  as  a  function  of  the 
parameter  /I.  In  all  cases,  the  curves  are  meant  to  denote  the  minimum  obliquity  at  which  ricochet  can  possibly 
occur  as  a  function  of  A.  In  Fig.  12b,  the  graph  is  presented  as  a  function  of  which  is  a  parameter 

proportional  to  the  striking  velocity.  The  abrupt  step  between  the  p>{A^—\y^^  plastic  curve  and  the  elastic 
P  =  {A^—iy^^  curve  is  predicted  by  the  model.  According  to  the  model,  there  can  exist  a  “dead  zone”  of 
intermediary  obliquities  for  which  the  conditions  of  neither  elastic  nor  plastic  ricochet  is  satished,  resulting  in 
the  stepwise  appearance  seen,  for  example,  in  Fig.  9. 
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10.  Conclusions 

The  long-rod  ricochet  algorithm  posited  by  Segletes  [1,2]  has  been  revisited  in  search  of  efficiency 
improvements.  The  original  method  required  an  iterative  approach  to  solving  the  equations  of  ricochet.  One 
key  step  in  accomplishing  this  improvement  is  the  replacement  of  a  moment-of-momentum  governing 
equation  in  the  original  model  with  a  simpler  minimization  constraint  on  the  rebound  angle  of  the  ricocheting 
projectile. 

With  this  change,  the  ricochet  equations  are  now  solved  in  analytical,  non-iterative,  form.  For  the  case  of 
elastic  ricochet,  there  are  three  branches  of  ricochet  behavior  depending  upon  the  relative  magnitude  of  the 
target  and  projectile  strengths,  projectile  density  and  striking  velocity.  The  critical  angle  of  elastic  ricochet  is 
quantifiable  as  a  function  of  the  elastic  ricochet  branch,  so  that  a  rapid  assessment  may  be  made  of  whether  an 
elastic-ricochet  solution  presents  itself.  In  the  absence  of  an  elastic  ricochet  solution,  a  non-iterative  algorithm 
is  employed  to  examine  particular  special-case  plastic  solutions  so  as  to  ascertain  that  solution  which  satisfies 
all  the  kinematic  and  force-limit  constraints,  including  that  for  rebound  minimization. 

In  the  new  algorithm,  no  iterative  steps  are  required  to  obtain  the  ricochet  determination/solution.  As  such, 
the  current  algorithm  detailed  in  Table  1  offers  a  much  more  streamlined  solution  than  that  embodied  in  the 
original  efforts  [1,2].  Minor  differences  in  results  between  the  two  methods  are  attributable  to  the  replacement 
of  the  moment-of-momentum  governing  equation  in  the  original  work  with  the  minimized-projectile-rebound- 
angle  constraint  in  the  current  effort.  Nonetheless,  the  revised  model  retains  key  elements  of  the  original 
model,  including:  ricochet  by  way  of  plastic  hinge  formation;  rod  rebound  from  the  target’s  ricochet  surface; 
target  gouging  as  a  means  to  redirect  the  target  force’s  line  of  action  during  the  ricochet  event;  and,  of  course, 
linear-momentum  conservation. 

Analytical  expressions  are  also  offered  for  the  critical  ricochet  obliquity  as  a  function  of  the  rod  and  target 
material  descriptions  and  the  striking  velocity.  These  expressions  are  summarized  in  Table  2.  In  some 
engagement  domains,  the  equations  are  exact  solutions  to  the  ricochet  equations,  while  in  other  domains,  the 
analytical  expressions  are  very  close  approximations  to  the  actual  solution,  acquired  through  the  use  of 
simplifying  assumptions. 


Appendix  A.  Derivation  aid 


The  math  presented  in  the  article  is  solely  based  on  solving  the  following  system  of  equations: 
Equation  in  A: 

cos(6  —  a)  cos  6 

a  — - , 

sm(d  -t-  q)  cos  q 

Equation  in  fl: 

_  cos(o!  -h  V]) 
sin(0  +  q) 

subject  to  the  material  constraints 

0<a^A,  q  —  0  (elastic  target) 
a  —  A,  qjt=0  (plastic  target) 

1  ^  /^eff  ^ 

the  kinematic  constraints  embodied  in  the  relation 


(A.l) 


(A.2) 


min[0,  njl  —  6  —  2q], 

as  well  as  the  a-minimization  constraint.  The  a-minimization  constraint  can  be  shown,  for  a  given  A  and  6,  to 
be  equivalent  to  the  solution  that  maximizes  /ieff. 

Derivation  of  key  equations  appearing  in  the  article  is  herein  described.  One  solution  technique  frequently 
used  in  this  appendix  involves  the  conversion  of  trigonometric  terms  to  double  angles,  to  eliminate  sine-cosine 
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products.  The  identities 

cos^  f;  =  (1  4-  cos  2q)/2 

(A.3) 

and 

sin  q  cos  q  —  sin  2q/2 

(A.4) 

are  employed  to  produce  equations  of  the  form  rco?,lr\-\-psmlr\-\- q  =  Q.  Isolating  the  cosine,  squaring  the 
equation,  and  converting  the  resulting  co^q  term  to  1— sin^j^  produces  a  quadratic  form  in  sinrj,  allowing 
solution  as 


sin  2t]  = 


—pq  —  r^y  +  p^  —  q^ 
r2  4-/72 


(A.5) 


Note  that  only  the  minus  root  is  retained,  so  as  to  keep  r]<d  for  the  range  of  q  values  to  which  this  equation  is 
applied  in  this  article.  The  value  of  q  is  then  obtained  by  way  of  an  inverse  sine.  Specihc  special-case  solutions 
are  now  addressed. 


A.  1.  Elastic  ricochet  (q  =  0) 

A.  1.1.  Branch  I 

In  branch  I,  it  is  the  value  of  A  which  limits  elastic  ricochet.  Eq.  (A.l),  as  a^A  with  q  =  0,  gives 
^  =cos(0— a)/tan  6.  One  may  observe  that  positive  a  values  raise  the  value  of  A  for  a  given  6.  Thus,  for  a  given 
A,  a  minimum  9  is  obtained  only  when  a  =  0.  When  a  =  0,  the  governing  relation  becomes  A=cos^6/sin6. 
Converting  the  cosine-squared  term  to  an  equivalent  sine  squared  term  and  solving  the  resulting  quadratic 
provides  the  branch  I  solution,  Eq.  (14). 

A.  1.2.  Branch  III 

In  branch  III,  it  is  the  value  of  fi  which  limits  elastic  ricochet.  Eq.  (A. 2),  as  t^  with  q  =  0,  gives  fi  =  cos  a/ 

sin  6.  Eor  the  a  =  0  limit,  this  becomes  jl=  1/sind  and  produces  Eq.  (17).  If  a  is  permitted  to  be  positive  and 
non-zero,  the  elastic  6^  may  be  reduced  further  for  a  given  /i.  However,  the  upper  limit  on  a  is  6.  And  so,  if  a  is 
permitted  to  take  on  this  upper  limit,  the  governing  equation  becomes  /i=  1/tan  0,  and  Eq.  (18)  follows 
directly. 

A.  1.3.  Branch  II 

In  branch  II,  the  initial  limiting  factor  is  the  value  of  fi,  and  thus  the  a  =  0  solution,  Eq.  (20),  is  identical  to 
that  for  branch  IIL  However,  as  a  is  allowed  to  increase,  the  value  of  A  becomes  the  limiting  factor,  prior 
to  the  point  where  a  increases  to  a  value  equal  to  6.  Thus,  for  the  elastic  6^,  both  a^A  must  be  satished, 
as  well  as  P-  Eq.  (A. 2),  with  these  limits  and  q  =  0,  gives  cos  a  =  /i  sin  9,  which  leads  to 

sin  a  —  V 1  —  cos2  a  =  \J \  —  ff'  sin^  9.  By  expanding  cos(d— a)  in  Eq.  (A.l),  one  obtains 

A  —  cos  a  cos2  0/  sin  0  4-  sin  a  cos  0.  (A. 6) 

Substitute  the  expressions  for  cos  a  and  sina  into  Eq.  (A. 6).  Isolate  the  square  root,  and  square  the  expression. 
Reduce  the  /l^(sin^  0  4-  cos2  0)  expression  to  leaving  just  an  expression  for  cos^  0,  which  can  be  solved 
directly  for  Eq.  (21)  (since  9  —  91  when  a  =  A  and  /ieff=  /!)•  The  limiting  a  value,  Eq.  (22),  is  obtained  from  the 
earlier  expression,  cos  a  =  /i  sin  0. 

A. 2.  Plastic  ricochet  ( a  =  A ) 

A.2.1.  Solution  for  a  specified  fleff 

Take  Eq.  (A.l),  solve  for  cos(0— a),  and  set  aside.  In  Eq.  (A. 2),  express  fi  +  q)  as  [(9  +  q)—{9—a.)], 
and  expand  as  a  cosine  of  the  sum  of  angles.  Eliminate  cos(0— a)  from  the  resulting  expression  by 
substituting  Eq.  (A.l),  and  solve  for  sin(0— a).  Take  these  two  resulting  equations  [one  for  cos(0— a) 
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and  the  other  for  sin(0— a)],  square  each  of  them  and  add  together.  The  result  should  be  an  equation  in 
which  o£  is  eliminated: 


l^lff  -  1  +  ^ 


2  cos^  r] 
cos^  6 


-  2/^eff^  ^^COS  (9 +  1])  =  0. 


Expand  cos{6  +  t]).  Clear  the  denominator.  Employ  the  technique  of  Eqs.  (A.3)-(A.5),  so  as  to  obtain  sin2>7 
and  hence  rj.  Eq.  (23)  is  analogous  to  Eq.  (A. 5),  with  the  terms  of  Eqs.  (24)-(26)  corresponding  to  the 
coefficients  that  comprise  Eq.  (A. 5). 

Many  expressions  can  be  used  to  obtain  a  expressions,  given  /leff  and  ?7/;eff,  for  example  solving  either 
Eq.  (A.l)  or  Eq.  (A. 2)  directly  for  a.  Elowever,  these  expressions  solve  in  terms  of  an  inverse  cosine,  which 
cannot,  unaided,  detect  the  presence  of  negative  values  of  a  (which  are  not  validly  permitted).  Instead,  start 
with  the  expression  which  was  earlier  in  this  section  derived  for  sin(d— a)  and  take  the  inverse  sine  to  isolate 
and  obtain  Eq.  (27). 


A. 2. 2.  Solution  for  which  a.  =  0 

Start  with  Eq.  (A.l).  Set  a  =  0.  Expand  sin(0  +  j^).  Isolate  cos^OjA  on  one  side  of  the  equation.  Employ  the 
technique  of  Eqs.  (A.3)-(A.5)  to  obtain  sin2j^  as 


.,  23  .  /I  A  sin  9  cos^  9 

sm  2ti  —  —  cos  9  —  sm  9  cos  9  —  2  sm  9  cos  9\  --\ - . 

'A  A  A^ 


{A.l) 


A  solution  to  Eq.  (A. 7)  will  exist  only  if  the  radicand  is  nonnegative.  Consider  this  radicand  inequality  to 
resolve  this  question  and  solve  for  the  allowable  values  of  A  (only  positive  A  solutions  are  relevant).  The 
resulting  inequality  may  be  expressed  in  the  form  of  Eq.  (28). 

Return  to  Eq.  {A.l)  to  solve  for  t].  In  two  places,  use  an  Eq.  (A. 4)  substitution.  Eq.  (29)  for  r\  follows. 

Eq.  (30)  follows  directly  from  Eq.  (A. 2),  when  a  is  set  to  zero. 

To  prove  the  assertion  that  no  ricochet  solutions  exist  if  Eq.  (28)  is  not  satished,  a  similar  derivation  to 
Eq.  (28)  may  be  obtained,  except  that,  instead  of  setting  a  to  0  at  the  outset,  it  should  be  set  to  a  positive 
constant  kq,  for  which  an  alternate  solution  is  sought.  The  derivable  inequality  for  this  condition  proves  to  be 

/  COS  0  \ 

2  —  AI - - - I  ^ 2  sin  0  (if  0 < a  =  ao ^ d  solution  exists). 

\cos(9  -  ao)/ 


But  if  Eq.  (28)  is  not  satished,  then  2  sin  0^(2— A)  is  implied,  leading  to  the  requirement  that 


cos  9  \ 

cos(9  —  ao)  J 


^2- A 


(A.8) 


be  satished  for  an  O<a  =  ao^0  solution  to  exist  simultaneously  when  there  is  no  a  =  0  solution.  Because  A  is 
positive,  Eq.  (A.8)  leads  to 


cos  9^  cos(9  —  ao). 


Elowever,  for  0<ao^9<n/2,  this  equation  can  never  be  satished.  One  may  conclude  therefore,  if  Eq.  (28) 
cannot  be  satished,  that  no  ricochet  solution  in  the  range  0  <  a  ^  6  can  exist. 


A.  2. 3.  Solution  that  maximizes  a 

Solve  Eq.  (A.l)  for  cos(0— a),  and  call  this  term  k.  Expand  sm{9  +  ri).  Employ  the  technique  of  Eqs. 
(A.3)-(A.5)  to  solve  for  sin2?7: 


sin  2t]  = 


(2k/ A)  —  tan  9  —  tan  9^  1  +  (4k  tan  9/ A)  —  (4k^/A^) 
tan^  9  +  1 


A  solution  is  possible  if  the  radicand  is  nonnegative.  The  maximum  allowable  a  is  9,  and  so  the  maximum 
value  for  k  is  unity.  When  k=l,  enforcing  a  nonnegative  radicand  gives  Eq.  (31).  Eor  the  case  where  the  k=  1 
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solution  exists,  substitute  k  =  1  and  make  use  of  trigonometric  identities  like  tan^  9+1=  sec^  6  to  obtain  the  rj 
solution  of  Eq.  (32).  Eq.  (33)  is  obtained  directly  from  Eq.  (A. 2)  for  the  case  where  a  =  9. 

Eor  cases  where  the  /:=  1  solution  is  not  viable,  solve  Eq.  (A.l)  for  cos(6— a).  Take  d/drj  of  the  equation  and 
set  equal  to  zero,  so  as  to  find  the  value  of  rj  for  which  cos(0— a)  is  maximized  [i.e.,  where  (9— a)  is  minimized]. 
The  solution  to  this  case  is  rj  =  l/2sin~*(cos  0)  =  ti/4—  9/2,  given  as  Eq.  (34). 

Define  as  k  the  cosine  of  this  minimized  angle,  so  that  /:  =  cos(0— o£)mm,  and  substitute  into  Eq.  (A.l);  solve 
for  k.  Substitute  r]  =  n/A  —  9/2  so  as  to  eliminate  t],  that  k  may  be  expressed  solely  in  terms  of  A  and  9.  Note 
that  cosr\  sm{9  +  r])  may  thus  be  reduced  to  1/2-  (1  +sin0),  to  obtain  Eq.  (36). 

To  obtain  Eq.  (37),  first  note  that  ((x.  +  ri)  =  —{9—(x.)  +  {9  +  t]).  Using  this  and  the  definition  of  k,  upon 
substitution  into  Eq.  (A. 2),  gives 


p=^l-k^ 


k 

tan(0  + 1])  ’ 


(A.9) 


Conversely,  the  definition  of  k  into  Eq.  (A.l)  gives,  with  some  regrouping, 


k 

tan(0  +  rj) 


A 

cos  9 


cos  rj  cos(9  +  t]) 


Substituting  Eq.  (A.  10)  into  Eq.  (A.9),  to  eliminate  k/tcin{9  + 1]),  gives 


(A.  10) 


R  =  V I  —  k^  -\ — ^^cos  M  cos(0  +  n).  (A.  11) 

cos  9 

Knowing  t]  =  n/4—9/2  allows  one  to  determine  that  cos  f/ cos(0  +  >^)  =  (cos  6)/2.  This  result  reduces 
Eq.  (A.ll)  to  Eq.  (37). 


A. 2. 4.  Solution  for  which  a.  =  ri 

Start  with  Eq.  (A.l),  and  isolate  cos(0— a).  Expand  sine  and  cosine  terms  for  the  angle  sums.  Multiply  the 
equation  by  sec^t],  so  that  all  t]  terms  appear  as  tangents  and  secants.  Square  the  equation,  substitute 
(1+tan^j^)  for  sec^t],  and  expand  the  terms  so  that  all  t]  terms  appear  as  powers  of  tanj^  terms.  Move  all 
the  terms  to  one  side  of  the  equation,  and  group  by  powers  of  tanj^.  Divide  the  equation  by  sin^0,  which  is 
the  multiplier  on  the  tan"^  rj  term,  in  order  to  obtain  a  fourth  order  equation  with  tan"*  rj  as  the  leading  term. 
Finally,  factor  out  a  cos  9  from  the  denominator  of  the  tan  fy  multiplier,  and  use  a  double  angle  formula  in  9,  to 
obtain  Eq.  (38). 

To  obtain  the  simple,  large-^  limit  to  rj,  Eq.  (39),  start  with  Eq.  (A.l)  and  approximate  all  cosine  terms  as 
having  a  value  of  1.  Solve  for  t]. 

To  obtain  the  accurate  approximation,  Eq.  (40),  start  with  Eq.  (A.l)  and  substitute  rj*  for  a.  Put  all  terms 
involving  t]  on  one  side  of  the  equation.  Expand  the  sm{9  +  r])  term.  Use  the  technique  of  Eqs.  (A.3)-(A.5)  to 
solve  for  sin2f7,  where  r=l,  />=  1/tan 0,  and  q=l—E,  using  the  term  grouping  for  E  in  Eq.  (41). 
A  straightforward  (sec^  0=1+  tan^  0)  trigonometric  substitution  is  needed  in  the  denominator  in  order  to 
obtain  Eq.  (40). 

Eq.  (42)  is  obtained  from  Eq.  (A. 2),  by  setting  a  to  r\. 

A. 2. 5.  Solution  for  which  r]  =  0 

From  Eq.  (A.l),  substituting  r]  =  d,  one  may  directly  obtain  the  manuscript’s  expression  for  Eq.  (44), 
by  noting  that  cos(0— a)  =  ^  tan 0.  To  obtain  the  expression  for  Eq.  (43),  combine  Eqs.  (A.l)  and  (A. 2) 
(with  r]  =  0)  to  obtain 

A  cos  a  r  4 

^“cos(0-o£)cos  0’  ^  ^ 

By  evaluating  cos  a  as  cos[0— (0— a)],  expanding  it,  and  utilizing  the  fact  that  cos(0— a)  =  A  tan  0,  the  expression 
for  Eq.  (43)  may  be  directly  obtained  from  Eq.  (A.  12). 
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Fig.  Bl.  Empirical  fit  for  critical  ricochet  obliquity,  with  (1  as  a  parameter,  applicable  for  0e>38.33°.  Fitted  curves  are  depicted  as  solid 
lines.  Actual  general-model  behavior  is  depicted  with  dashed  lines. 


Appendix  B.  Mode  IV-x  discussion 

Mode  IV-x  ricochet,  as  depicted  in  Fig.  5e,  plays  a  role  in  the  critical  obliquity  only  within  the  small 
triangular  (A,9c)  region  delimited  by  the  vertices  (1  /\/2,  45°),  (0.9373,  38.06°),  and  (0.8047,  42.497°).  To  show 
a  vertex  of  the  Mode  IV-x  region  at  (l/\/2,  45°),  insert  this  (A,0c)  point  into  Eq.  (38)  to  obtain  the  four  roots 
of  tanrja  =  ^  as  (0,0,— 1,-1).  The  double  roots  at  t7  =  0  indicate  a  terminus  of  the  Mode  IV-x  region.  Similarly, 
when  the  (A, 6c)  values  associated  with  the  vertex  (0.9373,  38.06°)  are  entered  into  Eq.  (38),  a  double  root  at 
r]=  17.31°  (^eff=  1)  is  produced,  indicating  the  other  Mode  IV-x  terminus  for  large-/!  ricochet.  At  the  third 
vertex  of  the  region,  (0.8047,  42.497°),  there  is  one  root  at  a  =  r]  =  0  and  one  at  r]=  15.83°  (^eff=  !)• 

If  one  were  to  attempt  an  empirical  characterization  of  the  Mode  IV-x  solutions,  it  proves  convenient  to 
note  (from  Fig.  5e  for  example)  that  the  a  =  t]  solution  associated  with  the  Mode  IV-x  limit  condition  (where  A 
is  lowered  to  the  point  where  the  a  and  rj  curves  are  tangent  at  a  single  point)  is  diminished  in  magnitude  with 
respect  to  the  a  =  >7  solution  associated  with  the  /!  =  1  condition.  We  therefore  employ  a  form  for  the  hnitely 
diminished  Kc  that  mirrors  the  form  of  the  perturbation  analysis  used  in  Eqs.  (55)-(56),  which  lends  itself  to 
ready  characterization  of  the  Mode  IV-x  solution  curve  for  large  /!,  as 


etc  —  »7c  =  30°  —  6c/3  —  Aa  (when  38.06°^  6c  ^6^=0) 
leading  to 


(B.l) 


cos((46c/3)  —  30°  -t-  Aa)  cos  6c 
sin((26c/3)  -b  30°  -  Aa)  cos((6c/3)  -  30°  -b  Aa) 


(when  38.06°^  6c  ^6«=o), 


(B.2) 


where  6^  =  0  is  the  angle  where  this  fitted  form  for  Mode  IV-x  solutions  merges  with  the  elastic  solution,  in  this 
case  at  6j(  =  o  =  45°.  The  near  linearity  of  the  diminishment  of  a  that  is  noted  in  Fig.  5e  suggests  that  the  Aa 
function  may  be  empirically  characterized  with  the  linear  approximation 


Aa  Rs  2.161(6c  -  38.06°)  (when  38.06°^  6c  ^45°). 


(B.3) 


Eq.  (B.3)  guarantees  that  ac  take  on  the  proper  value  at  6c  =  38.06°,  associated  with  the  /!  =  1  solution,  and  that 
it  takes  on  a  value  of  ac  =  0  at  6c  =  45°. 
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While  the  result  was  not  shown  in  Fig.  1 1,  it  lays  virtually  atop  the  Mode  IV-x  solution  curve  for  large  p.  In 
this  case,  large  /i  are  defined  as  being  in  excess  of  Vl,  since  the  upper  limit  of  this  fit,  at  0^  =  45°,  A  —  1/V2, 
which  blends  into  the  a  =  0  elastic  solution,  has  a  value  of  /ieff  equal  to  Vl  at  that  point. 

In  fact,  one  may  use  this  same  technique  to  empirically  approximate,  not  only  the  Mode  IV-x  solutions,  but 
all  the  limiting  (a^  =  f/c)  plastic  solutions  that  lie  between  the  /?  =  1  solution  and  the  branch  I  elastic  solution 
(/1>a/2),  for  which  no  analytical  function  is  otherwise  offered.  One  adopts  the  form  of  Eqs.  (B.1)-(B.2), 
utilizing  the  Aa  diminishment  term,  except  that  0^  =  0  is  no  longer  fixed  at  a  value  of  45°.  Rather,  the  upper 
limit  of  the  fit,  at  0^  =  0,  defining  the  point  where  the  fit  merges  with  the  elastic  solution,  varies  with  /?  according 
to  elastic  ricochet  relation. 


0,=o  =  sin-‘(l/i3eff).  (B.4) 

The  Aa  function  must  be  empirically  characterized  to  force  a^  to  a  value  of  zero  at  the  point  where  0c  =  0o(  =  o- 
Linearly  fit  to  this  constraint,  the  value  of  Aa  is  generalized  to 


Aa  = 


0c  -  38.06° 
0,=o  -  38.06° 


(30°  -  0,=o/3). 


By  setting  /leff=min  [/1,V2],  the  Mode  IV-x  fit  is  recovered  for  ^>41.  At  the  lower  limit  of  /!=  1,  Eq.  (B.4) 
indicates  that  0o,  =  o  =  9O°;  and  so  Aa=0,  which  allows  the  exact  fl=l  solution  to  be  recovered,  as  well. 

Eig.  B1  shows  the  quality  of  this  fit  for  predicting  the  critical  ricochet  obliquity  for  intermediate  values  of  p. 
In  the  figure,  the  fitted  predictions  are  drawn  with  solid  lines,  while  the  model’s  true  behavior,  acquired 
through  a  systematic  application  of  the  general  solution,  is  shown  with  dashed  lines.  One  may  observe  that  this 
empirical  fit  for  the  critical  ricochet  angle  does  an  excellent  job  of  matching  the  general  model  results,  over 
the  complete  range  of  p.  Because  this  fit  is  empirical,  however,  it  was  not  incorporated  into  the  contents  of 
Table  2. 
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308  VOYAGER  WAY 
HUNTSVILLE  AL  35806 

1  COMPUTATIONAL  MECHS 
CONSULTANTS 
J  A  ZUKAS 
PO  BOX  11314 

BALTIMORE  MD  21239-0314 

1  CYPRESS  INTERNATIONAL 
A  CAPONECCHI 

1201  E  ABINGDON  DR 
ALEXANDRIA  VA  22314 

2  DE  TECHNOLOGIES  INC 
R  CICCARELLI 
WFLIS 

3620  HORIZON  DR 
KING  OF  PRUSSIA  PA  19406 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


3  DOW  CHEMICAL  INC 
ORDNANCE  SYS 
C  HANEY 
A  HART 
B  RAFANIELLO 
800  BLDG 

MIDLAND  MI  48667 

3  DYNASEN 
J  CHAREST 
M  CHAREST 
M  LILLY 
20  ARNOLD  PL 
GOLETACA  93 1 17 

1  R  J  EICHELBERGER 
409  W  CATHERINE  ST 
BEL  AIR  MD2I0I4-36I3 

1  ELORET  INST 

NASA  AMES  RSRCH  CTR 
D  W  BOGDANOFF  MS  230  2 
MOFFETT  FIELD  CA  94035 

1  EXPLOSIVE  TECHLGY 
M  L  KNAEBEL 
PO  BOX  KK 
FAIRFIELD  CA  94533 

1  GB  TECH  LOCKHEED 
J  LAUGHMAN 

2200  SPACE  PARK  STE  400 
HOUSTON  TX  77258 

2  GB  TECH  LOCKHEED 
L  BORREGO  C23C 

J  FALCON  JR  C23C 
2400  NASA  RD  I 
HOUSTON  TX  77058 

6  GDLS 

38500  MOUND  RD 
W  BURKE  MZ436  21  24 
G  CAMPBELL  MZ436  30  44 
D  DEBUSSCHER  MZ436  20  29 
JERIDONMZ436  2I  24 
W  HERMAN  MZ  435  01  24 
S  PENTESCU  MZ436  21  24 
STERLING  HTS  MI  48310-3200 


3  GD  OTS 

D  A  MATUSKA 
M  GUNGER 
J  OSBORN 

4565  COMMERCIAL  DR  #A 
NICEVILLE  FL  32578 

2  GD  OTS 
DBOEKA 
NOUYE 

2950  MERCED  ST 
STE  13 1 

SAN  LEANDRO  CA  94577 

1  GENERAL  RSRCH  CORP 
TMENNA 
PO  BOX  6770 

SANTA  BARBARA  CA  93160-6770 

4  INST  FOR  ADVANCED  TECHLGY 
S  J  BLESS 

J  CAZAMIAS 
J  DAVIS 
H  D  FAIR 

3925  W  BRAKER  LN 
AUSTIN  TX  78759-5316 

I  INTERNATIONAL  RSRCH  ASSOC 

D  L  ORPHAL 
4450  BLACK  AVE 
PLEASANTON  CA  94566 

I  R  JAMESON 
624  ROWE  DR 
ABERDEEN  MD2I00I 

I  KAMAN  SCIENCES  CORP 
D  L  JONES 

2560  HUNTINGTON  AVE  STE  200 
ALEXANDRIA  VA  22303 

I  LOCKHEED  MARTIN  ELEC  &  MIS 
G  W  BROOKS 

5600  SAND  LAKE  RD  MP  544 
ORLANDO  FL  32819-8907 

I  LOCKHEED  MARTIN 
MIS  &  SPACE 
W  R  EBERLE 
PO  BOX  070017 
HUNTSVILLE  AL  35807 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


3  LOCKHEED  MARTIN 
MIS  &  SPACE 

M  A  LEVIN  ORG  81  06  BLDG  598 
M  R  MCHENRY 
T  A  NGO  ORG  81  10  BLDG  157 
1 1 1  LOCKHEED  WAY 
SUNNYVALE  CA  94088 

4  LOCKHEED  MIS  &  SPACE  CO 
J  R  ANDERSON 

W  C  KNUDSON 
S  KUSUMI 0  81  1 1  BLDG  157 
J  PHILLIPS  0  54  50 
PO  BOX  3504 
SUNNYVALE  CA  94088 

1  LOCKHEED  MIS  &  SPACE  CO 

R  HOFFMAN 
SANTA  CRUZ  FACILITY 
EMPIRE  GRADE  RD 
SANTA  CRUZ  CA  95060 

1  MCDONNELL  DOUGLAS 
ASTRNTCS  CO 

B  L  COOPER 
5301  BOLSAAVE 
HUNTINGTON  BEACH  CA  92647 

2  NETWORK  COMPUTING  SVC  INC 
T  HOLMQUIST 

G  JOHNSON 

1200  WASHINGTON  AVE  S 
MINNEAPOLIS  MN  55415 

1  PHYSICAL  SCIENCES  INC 
P  NEBOLSINE 

20  NEW  ENGLAND  BUS  CTR 
ANDOVER  MA0I8I0 

1  RAYTHEON  ELECTRONIC  SYS 
R  LLOYD 
50  APPLE  HILL  DR 
TEWKSBURY  MA  01876 

1  ROCKWELL  INTERNATIONAL 
ROCKETDYNE  DIV 
H  LEIFER 

16557  PARK  LN  CIR 
LOS  ANGELES  CA  90049 

1  SAIC 

M  W  MCKAY 

10260  CAMPUS  POINT  DR 

SAN  DIEGO  C A  92 12 1 


1  SHOCK  TRANSIENTS  INC 
D  DAVISON 

BOX  5357 

HOPKINS  MN  55343 

2  SOUTHERN  RSRCH  INST 
L  A  DECKARD 

D  P  SEGERS 
PO  BOX  55305 

BIRMINGHAM  AL  35255-5305 

I  ZERNOW  TECHNICAL  SVCS  INC 

L  ZERNOW 

425  W  BONITA  AVE  STL  208 
SAN  DIMAS  CA  91773 


ABERDEEN  PROVING  GROUND 

67  DIR  USARL 

AMSRD  ARE  SL 
P  TANENBAUM 
AMSRD  ARE  SL  BB 
R  DIBELKA 
EHUNT 
J  ROBERTSON 
AMSRD  ARE  SL  BD 
R  GROTE 
LMOSS 
J  POLESNE 
AMSRD  ARE  SL  BM 
D  FARENWALD 
DBELY 
G  BRADLEY 
M  OMALLEY 
R  SAUCIER 
A  DIETRICH 
AMSRD  ARE  SL  BS 
M  BURDESHAW 
AMSRD  ARE  WM 
T  W  WRIGHT 
AMSRD  ARE  WM  BD 
A  ZIELINSKI 
AMSRD  ARE  WM  MB 
W  DEROSSET 
AMSRD  ARE  WM  MD 
ECHIN 
G  GAZONAS 
J  LASALVIA 
AMSRD  ARE  WM  TA 
S  SCHOENFELD 
M  BURKINS 
N  GNIAZDOWSKI 
W  A  GOOCH 


NO.  OF 

COPIES  ORGANIZATION 


C  HOPPEL 
E  HORWATH 
D  KLEPONIS 
B  LEAVY 
MLOVE 
J  RUNYEON 
AMSRD  ARE  WM  TB 
R  BITTING 
J  STARKENBERG 
AMSRD  ARE  WM  TC 
R  COATES 
R  ANDERSON 
JBARB 
N  BRUCHEY 
T  EHLERS 
T  FARRAND 
M  FERMEN-COKER 
E  KENNEDY 
K  KIMSEY 
L  MAGNESS 
B  PETERSON 
D  SCHEFFLER 
S  SCHRAML 
A  SIEGFRIED 
B  SORENSEN 
R  SUMMERS 
W  WALTERS 
C  WILLIAMS 
G  RANDERS-PEHRSON 
AMSRD  ARE  WM  TD 
T  W  BJERKE 
S  BILYK 
D  CASEM 
J  CLAYTON 
D  DANDEKAR 
M  GRINFELD 
Y I  HUANG 
KIYER 
BLOVE 

M  RAFTENBERG 
E  RAPACKI 
M  SCHEIDLER 
S  SEGLETES 
T  WEERASOORIYA 
H  W  MEYER 
AMSRD  ARE  WM  TE 
J  POWELL 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


2  AERONAUTICAL  &  MARITIME 
RSRCH  LAB 
S  CIMPOERU 
DPAUL 
PO  BOX  4331 
MELBOURNE  VIC  3001 
AUSTRALIA 

1  DSTO  AMRL 

WEAPONS  SYS  DIV 
N  BURMAN  (RLLWS) 

SALISBURY 

SOUTH  AUSTRALIA  5108 
AUSTRALIA 

1  ROYAL  MILITARY  ACADEMY 
G  DYCKMANS 
RENAISSANCELAAN  30 
1000  BRUSSELS 
BELGIUM 

1  CANADIAN  ARSENALS  LTD 
P  PELLETIER 

5  MONTEE  DES  ARSENAUX 
VILLIE  DE  GRADEUR  PQ  J5Z2 
CANADA 

1  DEFENCE  RSRCH  ESTAB  SUFFIELD 
D  MACKAY 

RALSTON  ALBERTA  TOJ  2NO 

RALSTON 

CANADA 

I  DEFENCE  RSRCH  ESTAB  SUFFIELD 
C  WEICKERT 
BOX  4000  MEDICINE  HAT 
ALBERTA  TIA  8K6 
CANADA 

I  DEFENCE  RSRCH  ESTAB 
VALCARTIER 
ARMAMENTS  DIV 
R  DELAGRAVE 
2459  PIE  XI  BLVD  N 
PO  BOX  8800 

CORCELETTE  QUEBEC  GOA  IRO 
CANADA 

1  CEA/CESTA 
A  GEILLE 

BOX2LEBARP33II4 

FRANCE 


5  CENTRE  D’ETUDES  DE  GRAMAT 

C  LOUPIAS 
P  OUTREBON 
J  CAGNOUX 
C  GALLIC 
J  TRANCHET 
GRAMAT  46500 
FRANCE 

I  DAT  ETBS  CETAM 
C  ALTMAYER 

ROUTE  DE  GUERRY  BOURGES 

I80I5 

FRANCE 

I  FRENCH  GERMAN  RSRCH  INST 
P-Y  CHANTERET 
CEDEX  12  RUE  DE  LTNDUSTRIE 
BP  301 

F6830I  SAINT-LOUIS 
FRANCE 

5  FRENCH  GERMAN  RSRCH  INST 
H-J  ERNST 
F  JAMET 
P  LEHMANN 
KHOOG 
H  F  LEHR 

CEDEX  5  5  RUE  DU  GENERAL 
CASSAGNOU 
SAINT  LOUIS  68301 
FRANCE 

I  CONDAT 
J  KIERMEIR 
MAXIMILIANSTR  28 
8069  SCHEYERN  FERNHAG 
GERMANY 

I  TDW 
MHELD 
POSTFACH  13  40 
D  86523  SCHROBENHAUSEN 
GERMANY 

I  DIEHL  GBMH  AND  CO 
M  SCHILDKNECHT 
FISCHBACHSTRASSE  16 
D  90552  ROETBENBACH  AD  PEGNITZ 
GERMANY 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


4  ERNST  MACH  INSTITUT 
V  HOHLER 
E  SCHMOLINSKE 
E  SCHNEIDER 
K  THOMA 
ECKERSTRASSE  4 
D-7800  FREIBURG  I  BR  791  4 
GERMANY 

3  FRAUNHOFER  INSTITUT  FUER 

KURZZEITDYNAMIK 
ERNST  MACH  INSTITUT 
H  ROTHENHAEUSLER 
H  SENF 

E  STRASSBURGER 
KLINGELBERG  I 
D79588  EFRINGEN-KIRCHEN 
GERMANY 

3  FRENCH  GERMAN  RSRCH  INST 
G  WEIHRAUCH 
R  HUNKLER 
E  WOLLMANN 
POSTFACH  1260 
WEIL  AM  RHEIN  D-79574 
GERMANY 

2  lABG 

M  BORRMANN 
H  G  DORSCH 
EINSTEINSTRASSE  20 
D  8012  OTTOBRUN  B  MUENCHEN 
GERMANY 

1  INGENIEURBUERO  DEISENROTH 
AUF  DE  HARDT  33  35 
D5204  LOHMAR  I 
GERMANY 

1  NORDMETALL  GMBH 
L  W  MEYER 
EIBENBERG 
EINSIEDLER  STR  18  H 
D-09235  BURKHARDTSDORF 
GERMANY 

2  TU  CHEMNITZ 

L  W  MEYER  (X2) 

FAKULTAET  FUER  MASCHINENBAU 

LEHRSTUHL  WERKSTOFFE  DBS 

MASCHINENBAUS 

D-09I07  CHEMNITZ 

GERMANY 


5  RAFAEL  BALLISTICS  CTR 
E  DEKEL 

Y  PARTOM 

G  ROSENBERG 
Z  ROSENBERG 

Y  YESHURUN 
PO  BOX  2250 
HAIFA  3 1 02 1 
ISRAEL 

I  TECHNION  INST  OF  TECH 
FACULTY  OF  MECH  ENGNG 
S  BODNER 
TECHNION  CITY 
HAIFA  32000 
ISRAEL 

1  ESTEC  CS 

D  CASWELL 
BOX  200  NOORDWIJK 
2200  AG 
NETHERLANDS 

2  PRINS  MAURITS  LAB 
H  J  REITSMA 

E  VAN  RIET 
TNO  BOX  45 
RIJSWIJK  2280AA 
NETHERLANDS 

I  TNO  DEFENSE,  SECURITY  AND 
SAFETY 
R  ISSELSTEIN 
PO  BOX  96864 
THE  HAGUE  2509JG 
NETHERLANDS 

I  ROYAL  NETHERLANDS  ARMY 
J  HOENEVELD 
VDBURCHLAAN3I 
PO  BOX  90822 
2509  LS  THE  HAGUE 
NETHERLANDS 

I  INST  OF  CHEMICAL  PHYSICS 

A  YU  DOLGOBORODOV 
KOSYGIN  ST  4  V  334 
MOSCOW 
RUSSIA 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


3  INST  OF  CHEMICAL  PHYSICS 

RUSSIAN  ACADEMY  OF  SCIENCES 
G  I  KANEL 
S  V  RAZORENOV 
A  V  UTKIN 

142432  CHERNOGOLOVKA 
MOSCOW  REGION 
RUSSIA 

1  INST  OF  EARTHS  CRUST 
P  I DOROGOKUPETS 
664033  IRKUTSK 
RUSSIA 

3  INST  OF  MECHL  ENGRG  PROBLEMS 

V  BULATOV 
D  INDEITSEV 

V  MESCHERYAKOV 
BOLSHOY,  61,  V.O. 

ST  PETERSBURG  199178 
RUSSIA 

1  INST  OF  MINEROLOGY  & 
PETROGRAPHY 

V  A  DREBUSHCHAK 
UNIVERSITETSKI  PROSPEKT,  3 
630090  NOVOSIBIRSK 
RUSSIA 

2  IOFFE  PHYSICO  TECHNICAL  INST 
DENSE  PLASMA  DYNAMICS  LAB 
E  M  DROBYSHEVSKI 

A  KOZHUSHKO 

ST  PETERSBURG  194021 

RUSSIA 

1  IPE  RAS 

A  A  BOGOMAZ 
DVORTSOVAIA  NAB  18 
ST  PETERSBURG 
RUSSIA 

2  LAVRENTYEV  INST 
HYDRODYNAMICS 
L  A  MERZHIEVSKY 

V  V  SILVESTROV 
630090  NOVOSIBIRSK 
RUSSIA 

1  MOSCOW  INST  OF  PHYSICS  &  TECH 
S  V  UTYUZHNIKOV 
DEPT  OF  COMPTNL  MATH 
DOLGOPRUDNY  1471700 
RUSSIA 


1  RSRCH  INST  OF  MECHS 

NIZHNIY  NOVGOROD  STATE  UNIV 
A  SADYRIN 

P  R  GAY  ARINA  23  KORP  6 
NIZHNIY  NOVGOROD  603600 
RUSSIA 

1  SAMARA  STATE  AEROSPACE  UNIV 
L  G  LUKASHEV 
SAMARA 
RUSSIA 

1  TOMSK  STATE  UNIV 
A  G  GERASIMOV 
5-THARMY  ST  29-61 
TOMSK  634024 
RUSSIA 

1  DYNAMEC  RSRCH  AB 
A  PERSSON 

PO  BOX  201 

SE-15123  SODERTALJE 

SWEDEN 

7  FOI 

SWEDISH  DEFENCE  RSRCH 
AGENCY 

GRINDSJON  RSRCH  CENTRE 

L  GUNNAR  OLSSON 

B  JANZON 

G  WIJK 

R  HOLMLIN 

C  LAMNEVIK 

LEAST 

M  JACOB 

SE- 14725  TUMBA 

SWEDEN 

2  SWEDISH  DEFENCE  RSRCH  ESTAB 
DIV  OF  MATERIALS 

S  J  SAVAGE 
J  ERIKSON 
S-17290  STOCKHOLM 
SWEDEN 

2  K&W  THUN 
WLANZ 
W  ODERMATT 
ALLMENDSSTRASSE  86 
CH-3602  THUN 
SWITZERLAND 


NO.  OF 

COPIES  ORGANIZATION 


5  DERA 
I  CULLIS 
IP  CURTIS  QI 3 
AHARTQI3 
KCOWANQI3 
M  FIRTH  R3I 
FORT  HALSTEAD 
SEVENOAKS  KENT  TNI4  7BP 
UNITED  KINGDOM 

1  UK  MINISTRY  OF  DEFENCE 
G  J  CAMBRAY 

CBDE  PORTON  DOWN  SALISBURY 
WILTSHIRE  SPR  OJQ 
UNITED  KINGDOM 

7  INST  FOR  PROBLEMS  IN 

MATERIALS  SCIENCE 
S  FIRSTOV 
B  GALANOV 
O  GRIGORIEV 

V  KARTUZOV 

V  KOVTUN 

V  MILMAN 

V  TREFILOV 

3  KRHYZHANOVSKY  STR 
252142  KIEV- 1 42 
UKRAINE 

1  INST  FOR  PROBLEMS  OF  STRENGTH 

G  STEPANOV 

TIMIRYAZEVSKAYU  STR  2 

252014  KIEV 

UKRAINE 


